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Abstract 

Homogeneous relaxation is a ubiquitous phenomenon in semiclassical kinetic theo- 
ries where the quasiparticles are distributed uniformly in space, and the equilibration 
involves only their velocity distribution. For such solutions, the hydrodynamic variables 
remain constant. We construct asymptotically AdS solutions of Einstein's gravity dual 
to such processes at strong coupling, perturbatively in the amplitude expansion, where 
the expansion parameter is the ratio of the amplitude of the non-hydrodynamic shear- 
stress tensor to the pressure. At each order, we sum over all time derivatives through 
exact recursion relations. We argue that the metric has a regular future horizon, or- 
der by order in the amplitude expansion, provided the shear-stress tensor follows an 
equation of motion. At the linear order, this equation of motion implies that the met- 
ric perturbations are composed of zero wavelength quasinormal modes. Our method 
allows us to calculate the non-linear corrections to this equation perturbatively in the 
amplitude expansion. We thus derive a special case of our previous conjecture on the 
regularity condition on the boundary stress tensor that endows the bulk metric with 
a regular future horizon, and also refine it further. We also propose a new outlook for 
heavy-ion phenomenology at RHIC and ALICE. 
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1 Introduction 

The AdS/CFT correspondence [JJ, fj, |3| has given us a framework to study non-equilibrium 
phenomena in gauge theories at strong 't Hooft coupling in real time. For conformal gauge 
theories at strong 't Hooft coupling, this correspondence, if applicable, implies that there 
exists a universal sector of non-equilibrium states. This universal sector of states maps 
to five-dimensional spacetimes with metrics which have regular future horizons, and are 
solutions of Einstein's equation with a negative cosmo logical constant. Construction of such 
solutions of gravity, perturbatively in the limit of slow spatial and temporal variations, leads 
us to uncover purely hydrodynamic phenomena in the dual gauge theory and also enables us 
to compute the hydrodynamic transport coefficients systematically ||, ||, ||, |7|, || || [10|, [11]. 
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It is, of course, of theoretical and possibly experimental interest, to uncover a wider class 
of non-equilibrium phenomena beyond hydrodynamics through the gauge/gravity correspon- 
dence. The universal sector of non-equilibrium states itself includes a huge spectrum of such 
non-equilibrium phenomena, including probably early time evolution of the quark-gluon 
plasma, as has recently been studied for boost invariant flows fl2"fl . [] 



The universal sector has a special characteristic in that all states constituting it can be 
uniquely characterized and their dynamics can be completely determined by the expectation 
value of the energy- momentum tensor alone ||14|| . This follows from the dual gravity de- 



scription. In an earlier work [15], we have proposed a field-theoretic explanation by drawing 
analogy to conservative solutions of the Boltzmann equation that we have constructed, and 
briefly explain below. 

It can be shown that an appropriate relativistic semiclassical Boltzmann equation cap- 
tures all perturbative non-equilibrium phenomena in non-Abelian gauge theories at suffi- 
ciently high temperature |16|, [Tj]]. The conservative solutions of the relativistic semicalssical 
Boltzmann equation exist for all values of the 't Hooft coupling and rank of the gauge group. 
These special solutions can again be completely characterized and their dynamics can be 
completely determined by the energy-momentum tensor alone. 

In these solutions, the components of the energy-momentum tensor follow a closed set of 
equations of motion which can be derived systematically from the Boltzmann equation. In 
addition to the equations for conservation of energy and momentum, this set also contains 
equations for the evolution of the shear-stress tensor. The other parameters of the quasi- 
particle distribution function, like the heat current for instance, do not decouple, but are 
algebraically determined by the energy-momentum tensor and their spatial derivatives in a 
local inertial frame where the mean velocity of the quasiparticles vanishes, and thus have no 
independent dynamical parts. 

The dynamics is determined by the energy-momentum tensor alone, because any solution 
to the closed set of equations of motion of the energy-momentum tensor can be lifted to 
a unique solution of the Boltzmann equation through the other algebraically determined 
parameters. 

A special class of solutions to these equations are purely hydrodynamic in nature and 
are known as normal solutions in the literature. f\ Further, any solution of the Boltzmann 



1 For related work with sources for field theory operators like metric perturbations turned on, please also 
see f|. 

2 The normal solutions were first found by Enskog [T^] in order to provide a systematic way of calculating 
transport coefficients from the Boltzmann equation. These solutions can be determined exactly by the 
hydrodynamic variables, as all other parameters of the quasiparticlc distribution are determined algebraically 
by the hydrodynamic variables and their spatial derivatives, in a local inertial frame where the mean velocity 
of the quasiparticles vanishes. Many formal aspects of these solutions were clarified in ||cj, and they 
were found for the relativistic semiclassical Boltzmann equation by Stewart plj . 
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equation can be approximated by an appropriate conservative solution at sufficiently late 
times [ p~5|l . 

We proposed that the conservative solutions, when extrapolated to strong 't Hooft cou- 
pling and large rank of the gauge group, exhibit universality in their dynamics and constitute 
the states of the universal sector. This explains why all states in the universal sector can be 
determined by the energy-momentum tensor alone. 

This proposal amounts to saying that the dual solutions of pure gravity will have regular 
future horizons, provided the boundary energy-momentum tensor follows a closed set of 
equations of motion. The constraints of the equations of motion of gravity imply conservation 
of energy and momentum. However these do not suffice to ensure regularity of the future 
horizon in the bulk, or determine the evolution of the boundary energy-momentum tensor 
completely. Our proposed set of equations also includes the equations of motion for the 
shear-stress tensor, which determine the evolution of the boundary energy-momentum tensor 
completely and ensure regularity of the future horizon in the bulk for the right values of the 
phenomenological parameters. 

The structure of these equations cannot be related to the corresponding structure in the 
case of conservative solutions because we do not have a known kinetic description at strong 
coupling. However, if we know the hydrodynamic transport coefficients up to some orders, 
we can phenomenologically construct the most general equation of motion of the shear-stress 
tensor in appropriate expansion parameters as we will discuss later. This equation has 
purely hydrodynamic solutions, just like conservative solutions, but it also has more general 
solutions corresponding to non-hydro dynamic relaxation. 

The aim of this work is to chalk out a course for derivation of equations of motion for the 
energy-momentum tensor, directly from the requirement of regularity of the future horizon 
in the gravity dual, for a special case of non-hydro dynamic relaxation. This special case is 
homogeneous relaxation. We will concentrate on the homogeneous case for two important 
reasons, firstly because they provide the simplest instances of purely non-hydro dynamic 
approach to equilibrium in general, and secondly because we can easily make a connection 
here with the physics of non-hydrodynamic branches of quasinormal modes of black branes. 

Homogeneous relaxation to equilibrium is a feature of special solutions of all kinetic 
theories. In the case of the Boltzmann equation, for instance , these solutions correspond to 
the quasiparticle distribution function being spatially uniform, but being non-Maxwellian in 
velocity space. The higher moments of the velocity distribution relax such that they vanish 
at equilibrium, while the first five velocity moments corresponding to the hydrodynamic 

3 Our study of homogeneous relaxation shows that the amplitude expansion of these equations of motion 
for the shear-stress tensor was not done correctly in It was actually misguided by the structure of 

the amplitude expansion in the Boltzmann limit. Consequently, we failed to connect our equation with the 
quasinormal modes after linearization. We will give the correct general phenomenological equation here. 
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variables - the density, the mean velocity, and the temperature - remain constant. Q 

In such solutions, the components of the energy-momentum tensor in the globally defined 
inertial frame with constant zero mean velocity, takes the form, 

^oo = e, hi = Uo = 0, Uj = p5ij + 7Ty(t), (1.1) 

where e is the energy density, p is the pressure dependent on the constant temperature 
T and 7Ty(t) is the shear-stress tensor. This form of the energy- momentum tensor trivially 
satisfies conservation of energy and momentum, dH^ = (or equivalently the hydrodynamic 
equations) for any 7Ty(i). We also require e = 3p and ir^Sy = so that the energy- momentum 
tensor is traceless. 

A special class of such homogeneous solutions are also conservative solutions, where the 
full solution can be determined from the solution to the closed set of equations of motion 
for 7f ij(t). These equations can be systematically expanded in the amplitude expansion 
parameter, which is the typical value of iiij divided by the constant pressure. 

Remarkably, the scalar channel of quasinormal modes for AdS black branes, which have 
zero wavelengths [26 , 27[ f28[ p9j , [30| [], also suggest a similar form of the energy-momentum 



tensor. [] We have a tower of such quasinormal modes, so that for such perturbations the 
energy- momentum tensor takes the following form by gauge/gravity duality, 

oo 

U*) = $ + E K^ e ~^ + <^ Rnt ) e-"**', (1-2) 

n=l 

where iff} corresponds to the unperturbed AdS black brane and ± UR n — iujj n (with uj n > 
0) are the non-hydrodynamic overtones in the scalar channel. Also a nfJiU is traceless and 
purely spatial. Given that t$ is traceless, it follows that a nfiu is traceless as well. As the 
perturbations are in the scalar channel, a nfJ , u is also transverse to the global boost four- vector 
m m of the black brane. It is not difficult to see that the form (|1.2Q of the energy-momentum 
tensor is equivalent to ( |1 . 1[ ) in the global inertial frame where the boost four- vector is (1, 0, 
0, 0). 

Here, we will develop a systematic method to construct solutions of gravity which have the 
boundary energy-momentum tensor of the form ( [Lip , by summing over all time derivatives, 



4 Spatially uniform solutions of the Boltzmann equation for special intermolecular potentials were first 
found by T. Carleman E2]. These solutions were also isotropic. More general solutions were later found by 



Wild pa] and Morgenstern |24 . 

5 For an earlier work please see p5[ . 

6 The quasinormal modes of AdS black branes reproduces all the poles in the retarded Green's function 
of the dual guage theory [ ^8| , , so they carry information regarding linear thermal relaxation in the gauge 
theory. 
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but perturbatively in the amplitude expansion parameter. [] 

We construct the solution both in Feffer man- Graham coordinates and the ingoing Eddington- 
Finkelstein coordinates. It will be easier to construct the solution in Feffer man- Graham 
coordinates first, because the constraints could be easily solved. The Feffer man- Graham 
coordinates will also be suitable for achieving a general construction of metrics, uniform in 
the spatial boundary coordinates. 

We can then obtain the metric in ingoing Eddington-Finkelstein coordinates by doing 
coordinate transformation systematically. The coordinate transformation also receives cor- 
rections order by order in the amplitude expansion, and at each order all time derivatives can 
be summed up. The whole method will be, in fact, a generalization of the method adopted 



in the purely hydrodynamic case in fl~4}| . 

The metric can be expected to be regular order by order in the amplitude expansion for 
appropriate shear-stress tensor 7Tij(t). The shear-stress tensor becomes arbitrarily small close 
to equilibrium, analogous to the spacetime variations of the hydrodynamic variables. In the 
latter case, the metric is regular order by order in the derivative expansion for right choice of 
transport coefficients. Thus we can expect that the metric should be regular for appropriate 
7Tij(t), order by order in the amplitude expansion after summing over all time derivatives at 
each order. As we will see, the time derivatives of ffij(t) will not become arbitrarily small 
even close to equilibrium. This will necessitate summing over all time derivatives at each 
order in the amplitude expansion. 

We can argue that the metric is manifestly regular in the ingoing Eddington-Finkelstein 
coordinates order by order in the amplitude expansion when 7Ty(i) follows an equation of 
motion. This equation of motion can also be found order by order in the amplitude expansion, 
summing over all time derivatives at each order. At the first order, we reproduce the scalar 
channel of quasinormal modes at zero wavelengths. At the next order, we obtain non-linear 
corrections, which can potentially modify the behavior drastically at early times. 

Organization of the paper : The organization of the rest of the paper is as follows. 
In Section 2, we will develop the method to construct solutions of gravity uniform in the 
spatial boundary coordinates, in the Fefferman- Graham system, order by order in the am- 
plitude expansion while summing over all time derivatives. In Section 3, we will translate 
these solutions to the Eddington-Finkelstein coordinates. In Section 4, we will analyze the 
regularity condition on the future horizon and see how it gives us an equation of motion 
for 7Ty(t). In Section 5, we will show how the analysis done here can be used to study the 



7 A class of translatioiially invariant and isotropic solutions were previously constructed in |3lJ, where 
the sources in the boundary corresponding to the non-normalizable modes were turned on and black hole 
formation in the bulk was studied. Our method of summing over time derivatives in the perturbation 
expansion here will be different and we can deal with arbitrary anisotropic configurations. We will adopt 
a different procedure for regularity analysis of the metric at the horizon here. For related earlier numerical 
study of homogeneous anisotropic asymptotically AdS metrics, please see [B2| . 
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space-time evolution of the matter created by ultrarelativistic heavy-ion collisions at RHIC 
and ALICE. Finally in the Discussion we will conclude by mentioning a few open questions. 



2 The metric in Feffer man- Graham coordinate system 
in the amplitude expansion 

In this section, we will develop a method to construct asymptotically AdS metrics, where 



the energy- momentum tensor at the boundary takes the same form (1.1) as in homogeneous 
relaxation in kinetic theories. The metric should have a regular future horizon when this 
boundary energy-momentum tensor satisfies further constraints, which will be investigated 
later. 

We will construct the metric order by order in the amplitude expansion, while summing 
over all time derivatives. We will use Fefferman- Graham coordinate system for both, a 
purpose of principle and a purpose of convenience. 

The purpose of principle is that the Fefferman-Graham coordinate system allows us to 
construct solutions for the most general form of the energy-momentum tensor corresponding 
to homogeneous relaxation, before we impose constraints for regularity at the future horizon. 



A similar point of view was taken in [[L4]] in the purely hydrodynamic case, where the metric 
in the bulk was constructed in Fefferman-Graham coordinates for purely hydrodynamic 
boundary energy-momentum tensor with arbitrary transport coefficients, in the derivative 
expansion. It was then shown that, for a unique choice of transport coefficients, the future 
horizon was regular order by order in the derivative expansion. 

In the present case of homogeneous relaxation, we have no standard phenomenological 
equation for the shear-stress tensor that can be reliably applied at strong coupling. Thus, 
in principle, we should construct the metric for the most general form of the shear-stress 
tensor at the boundary and then find out how the regularity condition at the future horizon 
in the bulk constrains it. When the boundary energy- momentum tensor is arbitrary, the 
Fefferman-Graham coordinate system is the suitable choice for constructing the metric. 

The purpose of convenience is that the Fefferman-Graham coordinate system allows us 
to satisfy the constraints easily. The constraints can be elegantly solved in the Fefferman- 
Graham coordinates, because they simply impose boundary conditions on the dynamical 
equations for evolution of the boundary metric in the radial direction. However, in other 
coordinate systems, the dynamical equations are not simply the evolution of the boundary 
metric. So, it is not easy to generate a simple algorithm to satisfy the constraints. We will 
find that the constraints here are more involved than in the purely hydrodynamic case, so 
the Fefferman-Graham coordinate system is indeed very convenient. 
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We will begin with obtaining the form of the metric order by order in the amplitude 
expansion. Then we will show how we obtain the equations for motion at each order in 
the amplitude expansion and how we can sum over all time derivatives at each order. The 
procedure will involve Fourier transformation of the time dependence, and this will introduce 
a few subtleties which we will also discuss. 



2.1 The form of the metric in the amplitude expansion 

Einstein's equation of gravity in five dimensions with a negative cosmological constant takes 
the following convenient form 

1 6 

Rmn — -RGmn = j^Gmn- (2.1) 

Any metric which solves the above equation will have constant scalar curvature R = — 4// 2 . 
Further, any asymptotically AdS metric which is a solution of this equation can be written 
in the Feffer man- Graham coordinate system in the form, 

ds 2 = j 2 (dp 2 + g^(p, z)dz»dz») . (2.2) 

This coordinate system should be able to cover a part of the the five- dimensional upper half 
plane p > 0, so it will have a coordinate or a real singularity only at a finite radial distance 



from the boundary p = [33, |H |. For the rest of this paper, we will choose our units such 
that 1 = 1. 

The boundary metric is defined as 

9{o)nv(z) = Yimg^ u (p,z). (2.3) 

By rules of gauge/gravity duality, which we will not review here, it turns out that the 
boundary metric is the metric on the four dimensional spacetime in which the gauge theory 
is living m. 

When the boundary metric is flat, i.e. when g^)^ (z) = r]^ u , the four dimensional metric 
9iw {Pi z ) i n ( |2.2j ) will have a Taylor expansion in the radial coordinate near p = |T3]] , such 
that 

oo 

g^(p, z) = r]^ + 9(2n)^u(z)p 2n . (2.4) 

n=0 

By the rules of gauge/gravity duality, when the boundary metric is flat, it also turns out that 
<7(4)/iK is the expectation value of the energy-momentum tensor in the dual non-equilibrium 
state (aside from a factor which is essentially a power of the rank of the gauge group |34], |36f| ) . 



7 



This result can also be obtained in a manifestly coordinate-independent manner [55, EE . 
Since we are in the limit where the rank of the gauge group is infinite, we will normalize 
tpv, the expectation value of the energy-momentum tensor in the dual state, such that 

We note that this identification of g^)^ wth t^ v makes sense only when the bulk metric 
has a regular future horizon. More generally, we should think of it as a boundary stress 
tensor. 

Einstein's equation automatically guarantees that g^)^ satisfies the equation of conserva- 
tion of energy and momentum, d^g^^ = 0, and is also traceless, i.e. it satisfies Tr(g^) = 



3^ , p6d - These allow us to identify it more generally with a conformally covariant energy- 
momentum tensor even when the bulk metric is not endowed with a regular future horizon. 

Einstein's equation for the metric ( [2.2|) reduces to a tensor, a vector and a scalar equation 
of motion. The tensor equation, which gives the dynamical equation of motion for evolution 
of the boundary metric in the radial coordinate, is 

\g" - |V - + -Tr{g- l g')g' - ±-Tr{g- l g>)g - Ric(g) = 0, (2.5) 

where ' denotes differentiation with respect to the radial coordinate p. The vector and tensor 
equations give the constraints on g^)^. The vector equation is 

V M TrGTV) - V%, = 0, (2.6) 
where V denotes the covariant derivative constructed from g. The scalar equation is, 

Tr[g- l g"] - ~Tr[<rV] - \t r\g^ g< g^ g>\ = 0. (2.7) 
Actually, it turns out that the true dynamical equation is 

\g" - |V - \g'g- l g' + -Tr^gy - Ric( g ) 



Q R(g) + -Tr(rV^V) - ^ (Tr(^V)) J = 0, (2.1 



+9 

which can be obtained by combining the tensor equation with the salar equation, after 



multiplying the latter with the four- dimensional metric g [0. Also, the true scalar constraint 
is 

R(g) + ?TrO/"V) + -Tr{g- l g'g- l g') - ^Tr^g')] 2 = 0, (2.9) 

which can be obtained by combining the trace of the tensor equation with the scalar equation. 
The best way to see how these equations work is to substitute the form of the metric 
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in (|2.S|) , fl2.6| ) and (|2.9p, and Taylor expand these equations in p about p = (the 
boundary). At the leading order, the dynamical equation fl2.8|) does not determine gu)^, 
but the vector constraint (|2.6| ) imposes the conservation of energy and momentum and the 
scalar constraint imposes the condition that it should be traceless. This is expected because 
5(4)^5 being the boundary stress-tensor, is indeed independent data other than the boundary 
metric gto)^ on the initial hypersurface p = 0, so it cannot be determined by the dynamical 
equation. On the other hand, the constraints impose the desired conditions on this initial 
data. 

Using the true dynamical equation fl2.8|) , we can determine all the higher coefficients in 
the Taylor expansion g^n) for n > 3 uniquely and algebraically in terms of polynomials of 
(7(4)^ and its derivatives, or equivalently in terms of t^ v and its derivatives. A few coefficients 
are 

1 

3(8),, = ^yv-^w^w)+^ n V - > (2-io) 

where all indices have been lowered or raised using the boundary metric rj or its inverse. It 
is easy to see why in these coefficients of the Taylor expansion (|2.4j) , the derivatives come 
only in pairs. The only sources of these derivatives in the dynamical eq. Q2.8Q are Ric(g) 
and R{g) where they do occur in pairs only. 

When the coefficients of the Taylor expansion (|2.4|) (like the ones above in (|2.10|) ), as 
determined by the true dynamical eq. ( j2.8|) , are substituted in the vector constraint fl2.6|) and 
scalar constraint ( j2.9| ), we find that the constraints trivially vanish, provided dH^ = and 
Tr(t) = (as obtained at the lowest order). The consistency of the power series expansion 
can be demonstrated once this trivialisation of the constraint equations has been proved to 
all orders and this has been done in M . 



The discussion above has been general and will apply to any t^ u as an initial data on the 
hypersurface p = 0. We will now specialize to the form of the energy- momentum tensor in 
a non-equilibrium state undergoing homogeneous relaxation to equilibrium. 

To specify this form, we have to define the hydrodynamic variables, the four velocity w M 
and the temperature T first. In kinetic theories, usually is defined as the local velocity 
of particle transport. However, it will convenient for us to use the so-called Landau frame 
in which is the local velocity of energy transport. By this definition, u M is also a timelike 
unit vector satisfying u^u^ = — 1. We will define the temperature such that the local energy 
density given by e = t^ v u^u v is (3/4)(7rT) 4 . We will also define b through b = \/{tiT). 

At equilibrium, of course, the hydrodynamic variables w M and T are constants, and the 
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energy-momentum tensor takes the following Lorentz-covariant form 



, _ "jU^Uy j- P^y M 111 



where is the tensor which projects to the spatial hyperplane orthogonal to w M , and is 
given by 

P flv = u fl u v + r] tlv . (2.12) 



It is easy to see that this energy-momentum tensor ( 2.11 ) is traceless. The equilibrium 
pressure is given by l/(46 4 ). 

The energy-momentum tensor of any non-equiibrium state can be parametrised by 

VM = 3 "' M ^; f - (z) + vM, (2.13) 

where the hydrodynamic variables are now functions of the boundary space and time coor- 
dinates, while 7r 'fw{z) is the shear-stress tensor which vanishes only at equilibrium. 

At non-equilibrium, the four-velocity u M becomes the local velocity of energy-transport 
and 3/(46 4 ) becomes the local energy-density. The shear-stress tensor tt^ should be such 
that it does not modify the local energy density and the energy current, so it must satisfy 
u^n^u = 0. Further the shear-stress tensor also satisfies Tr(ir) = 0, because the full energy- 
momentum tensor is traceless, and the equilibrium part is traceless by itself. Therefore, the 
shear-stress tensor ir^ actually has only five independent components, and together with 
the four hydrodynamic variables (the temperature T and three independent components of 
the four- velocity u M ), it can parametrize the nine independent components of a traceless 
boundary stress tensor. 

To specialize to homogeneous relaxation, we make the hydrodynamic variables constant 
and the shear-stress tensor dependent only on a global time, defined through —u^ = t, so 
that 

Given that u^n^ = 0, it is easy to check that dH^ = is satisfied for any ^^{t). 

We can further simplify the above form by going to the global comoving frame where 
= (1, 0, 0, 0). In this frame, we have, 

3 1 

^00 = toi — tiO — 0, tij = + 7r y(^)) (2.15) 



which is exactly the same as in ( [Lip with e = 3p = 3/(46 4 ). In this form, it is easy to see 
that the covariant tc^ indeed conforms with d^t^ = for any 7Ty(t). Also in this frame, 
the traceless-ness condition becomes ir^Sy = 0. Given that the conservation of energy 
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and momentum is trivially satisfied and also that the energy-momentum tensor is traceless 
for any traceless 7Ty(t), corresponding to any boundary stress tensor of the form (|2.14 ) or 
equivalently of form ( p.!5| ), a unique bulk metric is guaranteed to exist [14). Q 



The amplitude expansion parameter is the typical value of the shear-stress tensor tt^ u 
divided by the pressure l/(46 4 ). We note that the energy-momentum tensor which takes the 
form ( 2.14 ), or equivalently the form ( 2.15|) , is exactly given only by the sum of the zeroth 



order and the first order terms in the amplitude expansion, where the zeroth order term 
is the equilibrium part and the first order term is the shear-stress tensor itself. Unlike the 
purely hydodynamic case, the energy-momentum tensor itself does not get corrected order 
by order in the expansion parameter. 

At the zeroth order in the amplitude expansion, the metric is just that of the unperturbed 
boosted AdS black brane, dual to a finite temperature equilibrium state at zero chemical 
potentials in the gauge theory. The metric in Fefferman-Graham coordinates fl2.2| ) is fully 
specified by g^ v . We will refer to the AdS black brane g^ v as gj$, since it is the zeroth order 
term in the amplitude expansion. It is given by 



(4& 4 ~ p 4 ) 2 , A , p 4 

46 4 (46 4 + p 4 ) " u+ V 46 4 



C = - AZ,u4 r /^ u ^ + ( 1 + in ) ( 2 - 16 ) 



It is easy to see from the coefficient of the p 4 term that the boundary stress-tensor in this 
metric indeed takes the equilibrium form ( [2.11 ). The coefficients of the higher terms of 



the Taylor expansion about p = also agree with ( 2.10|) , where the derivatives should be 
dropped. 

In the global comoving frame, 



<® = - }£.Z?L , «ff=«S , = <>. ft?' = (l + £ K- (2-17) 



,(0) _ I 40 ~ P ) jo) _ JO) _ n „(0) 

46 4 (46 4 + p 4 )' 

At the first order in the amplitude expansion, the possible corrections to g^ v should 
involve tensors of rank two which have one tt^ in them. It is not hard to see that since w M is 
constant, and u^ir^ = by definition of n^ V) we get that u fJ, (u-d) n TT fJjU or u^idj 'dt) n 'K lJlv should 
vanish for all n. Further, because Tr(ir) = 0, the corrections to the metric at the first order 



in amplitude expansion can only be proportional to □ n 7T Mi ,, or equivalently (—) n (d/dt) 



2n 



TV,, 



Therefore, in the global comoving frame these corrections are of the form (—) n (d/dt) 2n 7Tij. 

An easier way to arrive at the conclusion above is to go to the global comoving frame 
directly, where we can readily see that the corrections take the form (d/dt) 2n 7Tij. We will 



8 It is not so obvious given that the boundary stress-tensor is not Cauchy data. But, it has been shown in 
jri) that a unique bulk solution exists if the bulk metric can be smoothly connected to the AdS black brane 
by turning on one or more perturbation parameters. 
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stick to the global comoving frame for the rest of this paper, until we extract the general phe- 
nomenological equation for the shear-stress tensor in Lorentz-covariant form. The regularity 
analysis of the metric is also convenient in the global comoving frame. 

Further inspection of the coefficients like fl2.10| ) of the Taylor expansion ( |2.4p of in the 
Fefferman-Graham metric ( |2.2| ) shows that, at the first order in the amplitude expansion, 
the corrections to g^ v should be of the form, 



oo oo , , x 2n 

n=0 m=0 X ' ' 



in the global comoving frame. Putting all our arguments together, at first order in the 
amplitude expansion, g^ v in the Fefferman-Graham metric ( |2.2|) takes the form, 

oo / i \2n 

9$ = <&> = 9% = 0, gff = £ b^f^(p) (| j (2.18) 

in the global comoving frame. 

By definition, the functions f^' 2n \p) are dimensionless, so they should be functions of 
the dimensionless variable {p/b). Also, 

/ / n \ 4+2n 



= ° ( ( £J ) > at p = °. ( 2 - 19 ) 



so that the vector constraint ( [2.6|) and the scalar constraint (|2.9| ) are satisfied. Explicit 
calculations of the coefficients like ( |2.10p , consistent with these constraints, clearly show 
that the leading term has to be coefficient of [dj d£) 2n ^iji or equivalently n n t^„. 



The dynamical equation ( p.8[ ) determines all the leading and subleading terms of jf( 1 > 2n ) 



for all n, except for the leading p 4 term of f^ 1 ' '. The correction to the boundary stress-tensor 
at the first order in the amplitude expansion is 7Ty(t), therefore its coefficient has to be p 4 . 
This implies 



/ a ' 0) = (£) 4 +o((£j )■«*<> = o- <^°> 



The subleading terms of Z^ 1 '^ are all of 0((p/b) in ). This is because the coefficients like ( p. 10 ) 



of Taylor expansion of g^ can pick up powers of 4n of p only in absence of derivatives, since 
p An accompany polynomials of the boundary stress tensor. A similar logic shows that /( 1,2n ) 
involves powers of 4 + 2n + 4m of p/b only, for m > 0. 

It will be convenient to do the Fourier transform of ^(t), so that, 



TTijit) = / du e'^iTijiuj). (2.21) 
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Similarly, we can do the Fourier transforms of the components of g$. Let us further define 
f^(p, to) such that 



(2.22) 



n=0 



With the above definition, g$ in Fourier space can be captured in terms of this single 
dimensionless function f^'(p, to). This single function encodes all the time derivatives of tt, 
since the Fourier transform of ( |2.18| ) implies 



The reality of 7Ty(t) and g$ implies that 

■ai-ui) = fa(u,))\ f {1) ( P ,-cu) = (/ (1) (p,^))*. 



(2.23) 



7T, 



(2.24) 



The latter is readily satisfied by the definition (|2.22 ) of f^(p,co), as the /^ ,2n ^(p)'s are real. 

The dynamical equation ( |2.8|) determines all coefficients of Taylor expansion of f^(p, oj) 
given that the coefficient of the leading term p 4 /6 4 is 1, so that the boundary stress tensor 
is corrected by n^. We thus find that 



p 4 b 2 u: 2 p 6 fl bW\ p 



12 W 



4 384 J b 8 



^ + ^)Q + ^). (2.25) 
30 23040/ b 10 KH ' V 1 



In the following subsection we will derive the equation of motion for f^(p,uj), in terms of 
which we capture the full metric at first order in the amplitude expansion. 

Similarly, we can argue that at second order in the amplitude expansion, the corrections 
t° SV * n ^ ne Fefferman- Graham metric ( |2.2|) should take the form 



(2) 
#00 



(2) 
(2) 

9ij 



oo n n / / 7 \ 2a / ? \ 26 



n=0 m=0 

= 0, 

oo n 



a,b=0 
a+b=n,\a— b\=m 



oo n n / / j \ 2a / , \ 2b 

EE^^w E (I) »■»(«) (I) 

n=0 m=0 " »— n \ V / V / 



a,fe=0 
a+6=ra,|a— b|=m 



+EE fo8+2n /i (Wm) ^) E 



n=0 m=0 



a,b=0 
a+b=n,\a— b\=m 



it) ^(i) ^ {t) 



i 6ij \ ,it 



1 „ ( d\ 2a , s ( & 



dt 



2b 



7T p? (£) 



.(2.26) 
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We can easily see that in the summation above over a and b for fixed n and m, when m ^ 0, 
(a, b) is either ((n + m)/2, (n — m)/2), or ((n — m)/2, (n + m)/2), and when m = 0, they 
are (n/2,n/2). We see that the dimensionless functions f^ 2 ' 2n ' 2rn ^ s capture corrections to 
the time-time components of the metric, f( 2 ' 2n ' 2m ^ s capture corrections proportional to the 
purely spatial identity matrix, and j| 2 ' 2n,2m )' s capture corrections which are purely spatial 
and traceless. 

In Fourier space, the above can be captured in terms of just three dimensionless functions 

/OO PCsC 
dbj e~ luJt I du x f!f\p,U),U 1 )'K p q(u 1 )'K pq (0U - Wi), 
-oo J —oo 

J2) _ (2)_Q 



OG 



lif = & % / du e luJt dui f { 2 ) (p,u,uj 1 )'K pq (uj 1 )'K m ( 

J —oo J —oo 

oo POO 

iut / j, , /(2), 



where 



2 

-g<%7r ra (wi)7r ra (u; - wi) , (2.27) 



71=0 m=0 

/f' 2n,2m) (p), for z = l, 2, 3. (2.28) 



It is important to note that the definitions of (p, w,wi) are such that they is symmetric 
under the exchange of cji and u; — wi. The powers of U\ above denote time derivatives acting 
on the first n and the powers of to — U\ denote time derivatives acting on the second tt, so it 
is natural to define these functions in the symmetric fashion mentioned. 

( 2) 

In the following subsections we will obtain equations of motion for these fl (p, to, Uxys. 
However, from the discussion so far, it is easy to see that the dynamical equation ( |2.8| ) 

(2) 

determines all the coefficients of f- (p, u, WiJ's in the Taylor expansion about p = 0, given 
the Taylor series for f^(p,tu) has been obtained previously at the first order. We thus 
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obtain, 



- ^ — ^ + ^IP J> 

(2 ) p 8 ^ 2 K 2 + (^-^i) 2 ) p 10 , 12 . 

/a " w~{ go 720 jva+OW )' 

/s OwO " 246^ " [ 240 240 + [p M } 



We can also check that the vector constraint (|2.6|) and the scalar constraint (|2.9|) vanish 
at the second order in the amplitude expansion, if we substitute the above expansions in 
them and Taylor expand about p = 0. We also note that the above expansions are 0(p 8 ) at 
the leading order, so that the boundary stress-tensor receives no corrections at this order as 
discussed before. 

(2) 

We note that the reality of g)J implies 

jfV-^-^i) = (jfV^i))*- (2-30) 



The above is guaranteed by the definitions ( [2. 28 



We can similarly obtain the forms of the expansion of g^ v at higher orders in the ampli- 
tude expansion. The corrections can be always grouped into three categories, namely the 
corrections in the time-time component, the corrections which are proportional to the spatial 
identity matrix, and finally the corrections which are purely spatial and traceless. In each 
category, we can have one or more independent tensor structures. However, after Fourier 
transforming the dependence on the time coordinate, all these corrections can be captured 
by a finitely few functions, which also efficiently sum over all time derivatives. 

2.2 The equations of motion and their solutions in the amplitude 
expansion 

We proceed to obtain the explicit equations of motion of the metric, order by order in the 
amplitude expansion, summing over all time derivatives. We will see that in the Fefferman- 
Graham coordinates, one can generate a simple algorithm to satisfy the constraints. 

In the previous section we made a distinction between the true dynamical equation (|2.8| ) . 
and the four-dimensional tensorial components (|2.5|) of Einstein's equation. The former 



was obtained by suitably combining the latter with the scalar component ( |2.7|) of Einstein's 



equation multiplied by the four-dimensional metric g. This distinction will not matter much 



in practice. The tensor equation ( |2.5| ) is almost as good as the true dynamical equation 
( [2.8| ) in determining the coefficients g^n)^ of the Taylor expansion (|2.4j) of g^ u for n > 2, 
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given a flat boundary metric <7(ow and a particular boundary stress-tensor which satisfies 
energy- momentum conservation, is traceless, and is identified with <?(4w. It is almost as 
good because it fails only in the case of one coefficient of the Taylor expansion, for n > 2, 
which is <?(8W- It will turn out that this will just introduce a slight modification in the 
general scheme, only at the second order in the amplitude expansion. 



We will use the tensor equation (2.S) therefore to obtain the equations of motion, order 
by order in the amplitude expansion Once we impose the correct boundary conditions for the 
equations of motion, the solutions will automatically satisfy the vector constraint ( |2.6|) and 
the scalar constraint (|2.9|) . Also, we will not make a distinction between checking the true 
scalar constraint ( |2.9| ) and the scalar component of Einstein's equation ( |2.7| ). The former is 
obtained by combining the latter with the trace of the tensor equation ( |2.5| ), which we are 
now treating as the dynamical equation. So, it is sufficient to check if the latter is satisfied. 

2.2.1 The metric at the first order in the amplitude expansion 

At the first order in the amplitude expansion, we recall that g^ u takes the form Q2.18|) . After 
Fourier transforming the dependence on the time coordinate, the metric takes the compact 
form ( 2.23 ), which involves a single function f^(p,co) defined in ( 2.22 ). We can readily find 



the equation of motion of f^(p, u) by expanding the tensor equation up to first order in the 
amplitude expansion and then Fourier transforming the time dependence. 
Let us define the differential operator through 



where, 



T (126 4 -p 4 )(46 4 + 3p 4 ) T 128p 6 6 4 
L ii(P) = 17714 -IwTU , ZT\ , L i2{p) - 



p(46 4 - p 4 ) (46 4 + p 4 ) ' ^ y, J (46 4 - p 4 ) (46 4 + p 4 ) 2 
4 (4& 4 + p 4 ; 
(46 4 - p 4 ) 2 



L M = ,L L - ( 2 -32) 



The equation of motion of f^(p,u) obtained from the Fourier transform of the tensor 
equation is ( |2~5| ) is 

D iM f (1) (j>,w) = 0. (2.33) 



At the first order in the amplitude expansion, the vector equation ( [2.6|) and the scalar 
equation (|2.7| ) identically vanish. It follows from the definitions of 7Tjj, that in Fourier space 
7Tij(uj)Sij = 0. A simple inspection reveals this makes the both the vector and scalar equations 
vanish. 
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However, in order that the solution can be extended at higher orders in the amplitude 
expansion, without conflicting with the constraints, and also because the correction to the 
boundary stress-tensor at first order in the amplitude expansion has to be 7Ty, we need the 



solution to the equation of motion ( |2.33| ) for f^'(p, u) to satisfy the boundary condition 



+ o 



, at p = 0. 



(2.34) 



One can readily see that the above boundary condition specifies the solution of f^(p,u) 
uniquely, by Taylor expanding eq. ( |2.33| ) in p about p = and checking this determines all 
the coefficients of p 2n of the Taylor expansion of f^(p,u) for n > 3. One can also check 
that this reproduces the known Taylor series of (|2.25|) of f^'(p,uj). 

To sum up, we have thus determined that at the first order the metric in Fourier space 
is given by f^'(p, u), which is the unique solution of eq. ( |2.33j ) with the boundary condition 
( |2.34j ). This solution sums up the time derivatives to all orders. 

Unfortunately, Q2.33j ) is not exactly solvable unless u = 0. We can, however, devise the 
following strategy. We use Q2.22j ) which relates f( l '(p,u) to the metric in real time. This 
relation actually is also the Taylor expansion of f^(p, u) in to about tu = 0. It is easy to see 
that ( |2.33D in conjunction with ( j2.22j ) implies that 



Dif 



(1,0) 



(2.35) 



and 



^ i/ (l,2n) = Mp) / (l, 2 (n-l)) ;forn > 1; 

b 2 



(2.36) 



where D is the differential operator Z} lw at tu = 0, i.e. 

d 2 . , , d 



D > = ^ + L "M^ + 2 -«- (2 - 37) 

with L n and L 12 as defined in ( |2.32j ). Also, the boundary condition Q2.34] ) implies Q2.19D 



and (ET2q ). 

The solutions of the homogeneous eq. ( |2.35|) are 



(2.38) 



The boundary condition ( p. 20 ) for implies that 



= -2s M = -2 ( 1 + Q .eg 



(2.39) 
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The eqs. ( |2.36p alongwith the boundary conditions imply the following recursion series 

f {1 > 2n) (p) = jk-si(p) f P dp w( S2iP ] ( ,, L 13 (p')f^Hp) 
b 2 Jo W(s l) s 2 ){p) 

I' d P w(* X{P \( >^ P)f {1 ' 2{n ~ l) \p))i for n > 1. ( 2 - 4 °) 

where W(s\, sz) is the Wronskian of si and S2. 

Through the exact recursion relations given by (|2.39| ) and (|2.40| ), we thus efficiently 



sum over all time derivatives in the metric at first order in the amplitude expansion, or 
equivalently define the coefficients of the Taylor series ( |2.22| ) of p l '{p,oj) about to — 0. 

2.2.2 The metric at the second order in the amplitude expansion 

At the second order in the amplitude expansion, we have seen that we can express in the 
Fefferman- Graham metric (|2.2|) in terms of three functions f^ 2 \p, oo, Co>i), for i = 1, 2, 3 after 
Fourier transforming the time dependence as shown earlier in (|2.27| ). 

The tensor equation ( |2.5| ) gives the equation of motion for all three /^'s. When we 
Fourier transform the dependence on time of this equation, to be consistent with the defini- 

(2) 

tions of fl s, the Fourier transform is done such that it is symmetric under the exchange 
of oj i and uo — U\ . 

The part of this equation of motion which is purely spatial and traceless at the second 

(2) 

order in the amplitude expansion, gives the equation of motion for /{ , 

D l0J fP (p, co, Wl ) = Sf (p, cj, wi) , (2.41) 
where Di w is defined through ( |2.31|) , and 



S{ 2) (p,^i) = R^(p)(J^f 1 \p,u 1 )^ 1 \p,u;-u 1 : 

+R ( ${p) (/ (1 Hp,^)^/ a Hp.w-^) + / (1, (p.w-wi)^/ (1) (p,wi) 

+ (^(p) + wi(w - Ux)R?l{p)) / (1) (p,^i)/ (1) (p^ - a*), (2.42) 



with 



2 ■ 



(46 4 + p 4 ) ' LZ vry (46 4 + p 4 ; 
(2) 64p 6 fr 4 (2) 166 s 

^13 (P) = (46 4 + p 4)3 ^ ^14 (P)- (4& 4_ p 4 )2 - ( 2 " 43 ) 
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(21 

Since / x captures the contributions to the metric which are purely spatial and traceless 
(just like the entire contribution of the metric at the first order in the amplitude expansion), 
it follows that its equation of motion should be given by the same differential operator 
acting on it, but with source terms on the right hand side. In fact, this structure of the 
equation of motion for the purely spatial and traceless parts is maintained to all orders in 
the amplitude expansion. 

The equations of motion couple f% and f 3 2 \ Let us define the following differential 
operators 

d 2 d 
D 2 iu = + L21 (p) q- + L 22 (p) + w% 3 (p) , 

d 

D22 = £24 (p) 7^ + £25 0), 

d 2 d 
Du = dp 2 + Lsi ^ ~dp~ + Ls2 ^ ' 
d 

Di2u = L 33 (p)— + L u (p)+u 2 L 35 (p), (2.44) 



with 



86 4 (126 4 + p 4 ) 4p 2 (486 8 + 8p 4 6 4 + 3p 8 



L 2 i(p) = --7JU IZY7II4 I J ' L 2 2 (p) 



p(46 4 -p 4 )(46 4 + p 4 )' (46 4 -p 4 )(46 4 + p 4 ) 2 ' 

46 4 (4fe 4 + p 4 ) (4fe 4 + p 4 ) 4p 2 (126 4 + p 4 

L 2s(p) = ~~ 7771 ITT - ' L 24(P) = -7777 77, ^25 (P) ~ 



(46 4 -p 4 ) 2 ' ^ yrj p(46 4 -p 4 )' ZJVr/ (46 4 -p 4 ) 2 ' 
_ 2(326 8 - 36p 4 6 4 - p 8 ) 4p 2 (1926 12 - 336p 4 6 8 - 60p 8 6 4 + p 12 ) 

-^31 (P) -7777 -J777I4 , „4\ ' L 32(P) 



£33 (P) 



p(46 4 -p 4 )(46 4 + p 4 ) ' ^ VF; (46 4 -p 4 ) 2 (46 4 + p 4 ) 2 

3(4fr 4 -p 4 )(16fe 8 + 24p 4 fe 4 + p 8 ) 
p(46 4 + p 4 ) 3 ' 



12p 2 (4fr 4 -p 4 )(166 8 + 24p 4 fe 4 + p 8 ) 12fe 4 
L 34(P) - (464 + p4) 4 > L35(p) -"(46 4 + p 4 )- (2 - 45) 



The part of the tensor equation ( |2.5| ) which is proportional to the purely spatial identity 
matrix gives 

AnoJ 2 (P> wi) + D22/ 3 (2) (P, w, wi) = S'f } (p, w, wi) , (2.46) 
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where 



(/ {1) (p,o; 1 )^/ (1) (p^-u; 1 ) + /( 1 )(p,u;-a; 1 )^/ (1) (P^i)) 
+ (4 2 3 } (p) + - u;a)i4J(p)) /«(p, i*)/^, w - (2.47) 



with 



3(46 4 + p 4 ) ' 22 yrj 3p(46 4 + p 4 ) 2 
(2) 16p 2 6 4 (12& 4 + p 4 ) (2) _ 166 8 

^23 (p) - 3 (4& 4 + p 4 )3 > ^4 (p) - 3 (46 4 -p 4 ) 2 - (2 - 48) 

The time-time component of the tensor equation ( |2.5|) gives 

L» 31 / 3 (2) (p, w, Wl ) + £W 2 (2) (P, w, wi) = Si 2) (p, w, Wl ) , (2.49) 

where 

4 2) (p,^i) = 4?(P) (/ (1> (p,^)^/ (1) (p,w-^) + / (1) (p,w-wi)^/ (1 Hp.^)) 

+ (4?(P) + ("? + (w - i*) 2 + ^(w - a*)) irfCp)) 
/«(p,a; 1 )/( 1 )(p,a;-a; 1 ), (2.50) 



with 

?(2) 



26 4 (46 4 - p 4 )(166 8 + 24p 4 6 4 + p 8 ) 
p(46 4 + p 4 ) 4 ' 



4? = _16p^(4 & 4 p 4 )(16^24p 4 & 4 + p 8 ) g) = °* (2.51) 



The eqs. (|2.41|) , ([2.461) , and ( |2.49|) sum the dynamical equations of motion. 



As discussed previously, to be consistent with the constraints, we require that 

/fV^i) = O ((£f) ,at p = 0. (2.52) 



The equation of motion ( p.41 ) determines the leading and the subleading terms of the ex- 
pansion. We require the same asymptotic behaviour of / 2 an d f% , however the eqs. ( ]2.46|) 



and (|2.49|) only determine the sum of the coefficients of p 8 of and f^ 2 \ We have already 
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discussed that the origin of this incapability is that the parent tensor eq. ( |2.5p from which 
these are derived is not the true dynamical equation given by ( p.8Q . 

This incapability is easily rectifiable, since the tensor eq. (|2.5|) is only incapable of 
determining the full coefficient of the p 8 term in the Taylor expansion of g^ u about p = 0, in 
terms of the boundary metric and boundary stress tensor. We should therefore use the true 
dynamical equation (2J3) to determine the coefficient of p 8 in Taylor series of either or 

(2) _____ _______ 

/ 3 , then the eqs. ( P-46 ) and ( |2.49 ) determine the other coefficient and all other subleading 



terms of both the expansions. So, the boundary conditions of /I and f^ z> are 



e(2) 



/ 2 (2) (p,a;,wi) 

/3 (2) (p,W,Wi) 



86 s 



O 



+ 
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at p = 0, 



at p = 0. 



(2.53) 



One can check that, with the boundary conditions (|2.52|) and (|2.53|) , the equations of motion 
( p.41| ), ( p. 46 ) and ( 2.49|) reproduce the known Taylor series (PD of /f^sthat are consistent 
with the constraints. Also, these boundary conditions make it explicit that the boundary 
stress-tensor is uncorrected at this order. 

We have also observed before that checking the constraints is equivalent to checking 
the vector eq. ( [2.6[ ) and the scalar eq. ( |2.7| ). We have reproduced these equations at 
the second order in the amplitude expansion in appendix A. These equations vanish if we 
substitute the solutions of the dynamical eqs. (|2.41| ), (|2.46|) and Q2.49|) , determined uniquely 
by the boundary conditions ( f2.52| ) and ( |2.53| ) in them. We can readily check this by Taylor 
expanding these equations in p about p = and substituting the Taylor expansions ( [2.29|) 
of fi 8 in them. 

To sum up, we have captured g^ u in the Fefferman-Graham metric ( |2.2| ) at the second 
order in the amplitude expansion (|2.27|) , through three functions //(/?, oj, for i = 1,2, 3, 
which are determined uniquely by the equations of motion (|2.41|) , (|2.46|) and ( p.49| ), along 
with the boundary conditions ( |2.52| ) and ( |2.53| ). This allows us to sum over all time deriva- 
tives at this order. We also note that all of these equations along with the constraints, after 
Fourier transforming the time dependence, are local in both the arguments u\ and u\, simply 
because these equations are local in time as well. 

Unfortunately again, this system of linear ODEs comprising ( |2.41|) , ( |2.46|) and ( [2.49| ), is 
not solvable exactly unless u — uj\ — 0. Therefore, we repeat our strategy at the first order 
to sum over the u and uj\ dependence through an exact recursion series. 

(2) 

We use the previously defined Taylor expansion of the f- (p, u,Uxys in oj and u)\ as 
defined in ( |2.28D , and use the full equations of motion (|2.41|) , ( [2.461 ) and ( p.49| ) to find 
equations of motion for the coefficients fi 2 ' 2n ' 2m \p) ■ We recall that, by definition < m < n. 
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The first equation of motion Q2.41 ) along with Q2.28| ) implies 

2 



(2,0,0) 



dp' 



dp' 



+R ( S(p) (/ (1 ' 0) ) 



(1,0)> 



(2.54) 



when n = m = 0, Di is as defined in (|2.37|) , 



(2,2n,2m) 



Ll3{P) ( f(2,2n-2,2m-2) 9 ,-(2,2n-2,2m) 



b 2 



+/i 



(2,2n-2,2m+2) 



0(n — m — 1) 
9(n - m - 2)) 



+_R (2) (p) A f(hn+m) j(l,n-m) 

11 <ip dp 



_|_^(2)/ \ | n(l,n+m) d j(l,n-m) _|_ y(l,n-m) j j(l,n+m) j 

12 \ (ip dp J 



+R ( S(p)f (1 ' n - m) f 



(l,n— m) x(l,n+m) 



n > 1, (m + n) = 0(mod 2), (2.55) 



and 



(2,2n,2m) 



^13 (p) 
6 2 



/ (2,2n-2, 2m -2) + 2 / (2,2n-2, 2m ) e(n _ m _ 1} 
+ / (2,2n-2, 2m+ 2) e(n _ m _ 2) j 

^(i,n-m-i)^(i,n+m-i) j n > i, ( m + n) = l(mod 2). (2.56) 



6 2 

The terms with Q(n — m — 1) contribute only when n > m + 1 and those with 0(n — m — 2) 
contribute only when n > m + 2. The R^s above have been defined in ( 2.43Q and L 13 has 
been defined in (gjg) . Also the /^'s are as obtained in fl2.39|) and fl2.40|) at the first order, 
and contribute only when p is even and non-negative, otherwise they vanish. 
The boundary condition (|2.52|) implies the solutions 



ft 



(l,2n,2m) 



(P) 



-si (jo) / dp 



S 2(p') 



o r W(s u s 2 )(p) 



/\ °1 



(2,2n,2m)/ ' 



+s 2 (p) / dp 



s i(p] 



o " W( Sl)S2 )(p') 



(2,2n,2m) 



(p), 



(2.57) 



where si and s 2 are the solutions of the homogeneous eq Dj/ = as defined in ( |2.38| ), 
S2) is their Wronskian, and gj 2,2 ™' 2 " 1 ) j s whatever that appears in the right hand 



side of equations of motion (|2.54|) , ( |2.55|) and ( |2.56| ) for /} ' n ' m ^'s, according to whether 
n = m = 0, or n > 1 and n + m = 0( mod 2), or n > 1 and n + m = 1( mod 2), 
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respectively. When n = m = 0, the integration above can be explicitly done, so that 



To find the equations of motion for f( 2 ' 2n ' 2m ^ s anc i ^ 2 > 2n > 2m )' s ^ ^ j g convenient to first de- 



couple the equations of motion ( p.46| ) and ( |2.49| ) for uj = uj\ = and then use the expansions 
(pi). 

Let us define the differential operator D 2 through 

~ dP_ (1445 s + 16p 4 6 4 + 3p 8 ) d 2 3(1 125 s + 16p 4 6 4 + 5p 8 ) d 
2 ~dp*~ p(46 4 + p 4 )(46 4 -p 4 ) d^ + p 2 (46 4 + p 4 )(46 4 -p 4 ) Tp ( ^ 

and the differential operator Z) 3 through 

~ d_ 4p 3 (12fe 4 + p 4 ) 
3 dp (46 4 + p 4 )(46 4 -p 4 )' 1 J 

The solutions of the homogeneous differential equation Z^/ = are 

The solution of the homogeneous differential equation D 3 f = is 

(46 4 - p 4 ) 2 . . 

u = ) A/ - " . 2.62 
46 4 (4fe 4 + p 4 ) v ; 

mi r r _ f (2,2n,2m) , , n(2,2n,2m) i ... 

ihe equations of motion tor j 2 s and j 3 can be written as 



,.(2,2n,2m) Q(2,2n,2m) 
2 — ^2 

.(2,2n,2m) c (2,2n,2m) 



^3/3 — *->3 5 (Z.bdJ 



where 



o( 2 , 2 n, 2m ) f d (48fe 8 + 8p 4 & 4 + 3p 8 ) \ (2) r ^ T ( 2 )^ 

52 - U~p(4& 4 +P 4 x4& 4 -p 4 )r 2 (p) ~ 24(p) 3 (p) ' 

^'"' 2m) = (^ (2) (p) " + ^^>°)^ + / 2 (2 ' 2 "' 2m) ) , (2-64) 
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4?(p) (|/ (1 ' 0) ) 2 + 24 2 2 ) (p) (/ (li0) |/ (1 ' 0) ) 

+#23 (p) (/ (1,0) ) 2 , when n = m = 0, 

^p) ^ / (2, 2 n-2,2m-2) + 2 / (2,2n^2, 2m ) e(n _ m _ 1} 

+/ 2 (2 ' 2n - W) 0(n-m-2)) 

+4?(p)^/ (1 " ) ^/ (1 ' n+m) 

+i? (2) (p) ( /( 1 ' n+m ) Ajf 1 ."-™) _|_ f (l,n-m) j(l,n+m) 

+i^?(p)/ (1 ' n_m) / (1 ' n+m) , wh en n > 1, (m+n) = 0(mod 2) 

LM ^ / (2,2n-2,2m-2) + 2 / (2,2n-2,2 m ) Q(n _ m _ 1} 

+/ 2 (2 ' 2n - 2 ' 2m+2) 0(n-m-2)) 

_£t 2 4^J /(lin _ m _ 1)/(liTt+m _ 1) ^ when = d 

b 2 



24? (p) (/ (1 ' 0) |/ (1 ' 0) 

+4 2 2 ) (P) (/ (1,0) ) 2 , whenn = m = 0, 

MP) ^ / (2,2n-2,2 m -2) + 2 / (2,2n-2,2 m ) @ ^ _ ^ _ ^ 
+ / (2,2n-2, 2m+ 2) @(n _ m _ 2) j 

+i? ( ' 2 ' ) (p) ( y( 1 ' n + m )_^_y( 1 ' n_m ) _)_ m) ^ j(l,n+m) j 

31 V rf P dp J 

+R < $(p)f {1 ' n ~ m) f {1 ' n+m) 

'(P) / y(l,n— m— 2) y(l,n+m) 

+ y(l,n-m)y(l,n+m-2)^ when ^ > ^ ( m+n ) = Q^ mod 3^ 
WP) ^ / (2,2n-2,2m-2) + 3 / (2,2n-2,2 m )@ (n _ m _ X ) 

+/ 2 (2 ' 2n - 2 ' 2m+2) e(n-m-2)) 



d(2)/ N 

• 3 ft 2 / (1 ' w " ro ~ 1) / (1,w+ro " 1) , when n > 1, (m+n) = l(mod 2). (2.66) 
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In the above eqs. i% s and R^'s are as defined in ( |2.48| ) and ( |2.48| ) respectively, and L23 
and L 35 are as defined in (|2.45| ). Further, the terms with Q(n — m — 1) contribute only when 
n > m + 1 and those with B(n — m — 2) contribute only when n > m + 2. The /( 1,p )'s are 
as obtained in ( |2.39| ) and fl2.40| ), and they vanish if p is not even and non-negative. 

The boundary condition ( |2.53|) determines the solutions of eq. ( [2.63| ) uniquely. The 



1 l- r r(2,2n,2rn) , 

solutions of J 2 s are 

, (2,0,0) ( , p 4 4fe 4 + p 4 46 4 + p 4 / ^46 4 - p 4 



66 4 126 4 V 4&4 + P 4 / 6&4 V \ 4&4 + P 4 
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a=0 



W a (p'^ 



= ^ t / dp / y ;7 v farn>l. (2.67) 



Here W(ti,t2,ts) denotes the Wronskian of t a 's defined in (|2.61|) . We recall that the 
Wronskian is the determinant of the 3x3 matrix whose i-th column is the transpose of 
(ti,dti/dp,d 2 ti/dp 2 ), for i = 1,2,3. Here W a denotes the determinant of the same matrix 
with the a-th column replaced by the transpose of (0,0, S^' 2n ' 2m ^), where gC 2 ' 2 ™' 2 ™) i s as 
defined in ( [2.64j ). 



The solutions of y^ 2 ' 2n ' 2m )' s are 

f(2 ,o,o) n (4& 4 -p 4 ) 2 2 (4& 4 -p 4 )(4& 4 + 3p 4 ) (4& 4 -p 4 ) 2 /46 4 -p 4 

h [P) 66 4 (46 4 + p 4 ) 3 (4fe 4 + p 4 ) 2 46 4 (46 4 + p 4 ) ° g V46 4 + p 4 

pp o(2,2n,2m) / '\ 

/i Wm) (p) = / dp' 5 * .,. (P) ,forn>l. (2.68) 
Jo u(p ) 



Here tt is as defined in (|2.62|) and gC 2 . 2 ™- 2 " 1 ) i s a s defined in ( [2.64j ). 



To summarize, the solutions ( p. 57] ), ( p.67|) and (|2.68| ) exactly capture the entire correction 



to the metric at second order, by specifying the coefficients of the Taylor series ( 2.28[ ) of 
fi\p,u,ui) in u and u)\ about u = u\ = 0, recursively. 

We can clearly extend our method to higher orders in the amplitude expansion, where 
we can solve the time dependence through similar recursion relations. 



2.3 Subtleties in the Fourier space 

We have seen that order by order in the amplitude expansion, all the time derivatives can 
be summed up efficiently in terms of a finitely few functions once we Fourier transform the 
time dependence. The equations of motion, which are ODEs in p with source terms, can 
also be obtained for these finitely few functions from the amplitude expansion of the tensor 
eq ( |2.5|) and Fourier transforming its time- dependence as well. 

We will point out here that the Fourier transform involves some subtleties. These sub- 
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tleties, are however, not new and could be readily illustrated with Navier-Stokes equation 
with no forcing term as an example. It is believed, though there is no rigorous proof, that 
the solutions will decay in the future and become homogeneous. Similarly, we believe in the 
case of AdS gravity, that the perturbations will decay and the solution will settle down to 
a boosted black brane. The point is that such decaying functions, like e~ at for instance, for 
a > 0, typically grow in the past and have no well-defined Fourier transform by the standard 
contour prescription, where the integration in t runs along the real axis from — oo to oo. 

The remedy for this is well known. It is convenient to state this remedy in terms of the 
inverse Fourier transform. Let 

= J- (_J^-Y a>0. (2.69) 
ztt \a — zoo J 

We see that f(oj) has a simple pole in the lower half plane at oj — —ia with residue —1 / (2iri). 
Let us inverse Fourier transform f(uj) by doing this integration over u using the following 
contour. We go slightly below the real axis and integrate in u from — oo — ie to oo + ie, 
where e is infinitesimal and positive, and then close the contour in the lower half plane by 
going along the semicircle at infinity. With this contour prescription, the inverse Fourier 
transform picks up contribution only from the pole in the lower half plane, so we get 

f(t) = e~ at , (2.70) 

which decays in the future. 

This contour prescription is simple, but has an important physical significance. We have 
a similar well-known contour prescription when one defines the retarded correlator in real 
time. The physical significance is that we pick up an arrow of time. In the case of solutions 
of gravity at the non-linear level, the result is even more drastic. 

We know that for purely hydrodynamic solutions of gravity, we can either have a regular 
future horizon or a regular past horizon, but not both jTJJ. The former happens at the first 



order in the derivative expansion when Ani]/ s = 1, and the latter happens when 4nr]/s = — 1. 
The regular perturbations of the AdS black brane and the AdS white brane cannot be 
glued without encountering some discontinuity, which is disallowed by Einstein's equation 
in vacuum with a negative cosmo logical constant. Therefore, the past of the solution which 
is purely hydrodynamic in the future will have drastically different behavior. 

The same situation will be repeated here in the case of homogeneous relaxation. In this 
case also, we will be interested in the solution which will settle to an AdS black brane in 
the future. Quite generally, the past behavior of the solution will depend on the specific 
initial condition. However, the behavior in the future can be expected to be governed by 
a generic phenomenological equation like the Navier-Stokes equation, provided the future 
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horizon is regular. Remarkably, this generic behavior at late time is exactly what is expected 
also in the dual theory. Our contour prescription isolates the generic future behavior from 
the non-generic behavior in the past, the latter being dependent on the specific details of 
the initial state. 

3 Translation to ingoing Eddington-Finkelstein coor- 
dinate system in the amplitude expansion 

The ingoing Eddington-Finkelstein coordinates are particularly suitable for analyzing the 
future horizon. The position of the future horizon at late time is revealed in these coordinates 
manifestly. In the case of stationary black holes, the metric is also manifestly regular at the 
future horizon in these coordinates. 

For more general geometries, which are perturbations of AdS black branes, we need to give 
a definition of the ingoing Eddington-Finkelstein coordinate system. This can be achieved 
by noting that the unperturbed boosted AdS black brane takes the following covariant form 



with P^ v being as in ( |2.12| ), the projection tensor on the spatial hyperplane orthogonal to 



We can think of r as the Eddington-Finkelstein radial coordinate and x^s as the boundary 
coordinates. The horizon can be located by finding where G^ v u^u u has a simple zero. We 
can readily see that the horizon is located at r = 1/6. Further the Hawking temperature 
T = l/(7r6), so we can legally use the same notation b as in the previous section, where b 
was defined as l/(7rT). The generators of the horizon in these coordinates are also limits 
of radial outgoing null rays. So this horizon is a future horizon. We also note that if we 
have chosen the outgoing Eddington-Finkelstein coordinates, the (fir) and (r/i) components 
of the metric would have been instead of — u^, and the horizon would have been a past 
horizon. 

The natural generalization of this form of the metric, to geometries which are perturba- 
tions of the boosted black brane, is 




(3.1) 



where is the boost parameter of the black brane as before, and 




(3.2) 




(3.3) 
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where w„ is now the local boost parameter and Gn V also depends on the boundary coordinates 



x M . As in ||14|| , we define the ingoing Eddington-Finkelstein coordinate system as the system 
in which the metric assumes the above form. 

When Vt? is a constant four-vector, as in the previous section, it is convenient to go 
to the global comoving frame where the boost parameter is (1,0,0,0). We define the 
ingoing Eddington-Finkelstein time v , as v — —u^ and call the spatial coordinates in the 
orthogonal spatial plane x*'s, where % = 1,2,3. The boosted black brane metric given by 
HI) and (13.21) now assumes the form 



ds 2 = 2dvdr + Gffi (r)dv 2 + Gf) (r)dzW, (3.4) 

where, 



Gtl = ~r 2 (l - , Gf = r%. (3.5) 

In our case of homogeneous relaxation, the metric in the global comoving frame in this 
coordinate system should assume the form, 

ds 2 = Idvdr + Goo( r ) v)dv 2 + Gij(r, v)dx l dx\ (3.6) 

where G 00 and G^ are also dependent on the time v. We note that G i0 = G 0i = as in the 
case of gio and goi in the Fefferman-Graham coordinate system, because there is no vector 
which can be constructed just from the boundary data 7Tjj(i). We will soon show that this is 
indeed the case when we obtain the metric explicitly by translating the Fefferman-Graham 
metric into this coordinate system. 

3.1 The equations for translation 

We have seen in the previous section that when the boundary stress-tensor assumes the 
form of homogeneous relaxation, the metric in Fefferman-Graham coordinates, in the global 
comoving frame, summing over all orders in the amplitude expansion, assumes the form 

ds 2 _ d P 2 + 9oo(p, t)dt 2 + gjjjp, t)dz l dz 3 

P 2 

If we do the following change of coordinates, 

P = *(r,v), (3.8) 
t = v + k(r,v), 
z l = x\ 
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the Fefferman-Graham metric ( |3.7p indeed translates to the ingoing Eddington-Finkelstein 
metric (|3.6|) provided 



-fr) +g 00 (p(r,v),t(r,v))^—J = 0, 

^ + «»Wr,«), t(r,,)) ^(l + ^J = * 2 , (3.9) 

where g o is now evaluated as a function of r and v. These equations thus determine 
the change of coordinates (|3.8|) from the Fefferman-Graham coordinates to the Eddington- 



Finkelstein coordinates. 

Once the Fefferman-Graham coordinates are determined as functions of the ingoing 
Eddington-Finkelstein coordinates, one can readily translate the Fefferman-Graham metric 
to the ingoing Eddington-Finkelstein metric. We note that, just like the Fefferman-Graham 
metric (|3.7|) is fully specified by g 00 and gij, the ingoing Eddington-Finkelstein metric (|3.6| ) 
is also fully specified by Goo and G^. The change of coordinates (|3.8[), determined by (|3.9|), 
implies 

d ®Y / / \ , / J, dk\ 2 „ _ 2 
-t^) +9oa(p(r,v),t(r,v))ll + - = G 00 $ 2 

gij (p(r,v),t(r,v)) = G %] e. (3.10) 

The above equations thus translate the metric. We note that the equations above are alge- 
braic equations and not PDEs like ( |3.9|) , where r and v derivatives of the unknowns $ and 
k appear. In the above equations, the r and v derivatives of the unknowns, Goo and do 
not appear. 

For the purely hydrodynamic case too, similar equations like ( |3.9|) and ( p.lOj ) for changing 
coordinates and translating the metric respectively from the Fefferman-Graham to ingoing 



Eddington-Finkelstein system were derived [|14j]. However, there the equations for changing 
coordinates and that for translating the metric did not disentangle nicely like here. 

For the unpertubed AdS black brane, &(°>(r, v) and k^(r,v), which completely specify 
the change of coodinates ( |3.8| ) translating the g^ u in the Fefferman-Graham metric ( [2.17D to 
the G^ v in Eddington-Finkelstein coordinates (|3.5|) are, 



$(°)(r,«) = $ {0) (r: v2/; 



v / & 2 r 2 + v ^4 r 4 _ I ' 

k {0 \r,v) = k {0) (r) = b(n-2 arctan(6r) + ^ log f J J . (3.11) 
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We can readily see that these solve the eqs. (PT5|) . Further, the eqs. ( |3.10D are also satisfied 
when G^q and G^ are given by ( |3.5| ). We also observe that both $(°)( r ) and k^(r) behave 
like 1/r, as we go towards the boundary at r = oo. 

We will now see how to implement eqs. ( |3.9| ) and ( |3. 1U[ ) for changing the coordinates 



and translating the metric respectively, from the Feffer man- Graham system to the ingoing 
Eddington-Finkelstein system, order by order in the amplitude expansion. At each order, 
we will sum over all derivatives with respect to the Eddington-Finkelstein time v through 
recursion relations, just like we summed over all derivatives in Fefferman-Graham time t in 
the previous section. 

3.2 Translation at the first order in the amplitude expansion 



At each order in the amplitude expansion, it will be easier to solve the eqs. ( |5T5| ) specifying 
change of coordinates ( |3.8| ) first. At the first order in the amplitude expansion, the change 
of coordinates do not receive any correction. We note that the change of coordinates from 
the Fefferman-Graham system to the ingoing Eddington-Finkelstein system, given by (3~£) 



maintains spatial translational and rotational symmetry of the former manifestly, to all 
orders in the amplitude expansion. It is easy to see that since itijSy = 0, there is no scope to 
write anything which is invariant under spatial rotations and translations, that can possibly 
correct $(r, v) and k(r,v), from their forms at the zeroth order given by ( |3.11| ). 



We recall that g^ v in the Fefferman-Graham coordinates at the first order in the derivative 
expansion takes the form ( ^.18[ ), where vanishes and g^' is traceless. The equations for 
translation of the metric ( |3.10| ) at the first order in the amplitude expansion therefore implies 



G 00 - 



(i) 

(1) _ gff (p = & °)(r),t = v + kW(r)) 
($(°)(r)) 5 



= yii KH ~ z: ~r (3.i2) 



where &°'(r) and k^(r) are as in ( |3.11| ) for the unperturbed AdS black brane. 



Further, just like g$ in ( |2.18|) , we can write as follows 



oo / , v n 

G^ = J2b 2+n F^\r)lf) 7r». (3.13) 

n=0 ^ ' 

The difference from ( p,18| ) here is that TTi j IS cL function of the Eddington-Finkelstein time 
v and not a function of the Fefferman-Graham time t. Additionally, the time derivatives 
run over both even and odd powers, unlike even orders only in the latter case. Also, by 
definition, F^ 1 ' n ' , s are dimensionless functions like f^ 1,2n 's. 
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Substituting flngp and ( g^ ) in ( gTg ), we get 



_ 

F(1 ' n) W = 7^^E/ (1 ' 2m Up = ^ (0 Hr))^ ( ° )(r)) 



n-2m 



($(°)(r)) ^ b n - 2m (n - 2m)! 

when n is even, 

n-l 

^ 2 S^i^W^Co^ (* (0) M) 



n—lra 



F (l.n)( r) = ? _^ / (l,2 m) ^ = $ (0) (r) ) 



($(°)(r)) ^ b n - 2m (n - 2m)! 

when n is even. (3-14) 

We recall that we have solved /( 1,2m ) in the previous section recursively and our solutions 
are given by (|2.39| ) and ( |2.40| ). Thus the equations above also give us recursive solutions for 

i* 1 ' 1 ' 71 ) (r) . We will not make these recursion relations more explicit here. The leading term 

F (i,o) 

(r) in this recursion can be easily found to take the form 

iXi.o)( r ) = _ftV log(l--LY (3.15) 



The recursion series (|3.14|) allows us to sum over all the time derivatives at the first order in 



the amplitude expansion. 

When we derived the recursion series in the Feffer man- Graham coordinate system to sum 
over all time derivatives, we went to Fourier space. This has not been necessary at the first 



order in the amplitude expansion here, because we had to only solve the algebraic eq. ( |3.12| ). 
At the second order in the amplitude expansion, the equations for changing coordinates will 
be PDEs and not algebraic equations, so Fourier transforming the time dependence will be 
necessary again. 

Therefore, for the sake of completeness, let us see if the metric here also assumes a 
compact form here in Fourier space too. The Fourier transform of itij{v) can be defined 
through 

/oo 
dco e-^itijlu). (3.16) 
-oo 

We note that since ^ij{v) is the same function of v as TTij(t) is a function of t. Given that we 
defined the Fourier transform of 7Ty(t) in ( |2.21| ) in the same way, it follows that the Fourier 
transforms TTij(uj) are identical. 

Similarly, we can define the Fourier transform of G^\r,v) in v. Let us also define 
fW(r,w) in analogy with f^(p,u) in Q2.22| ) ,as 



FW(r,w) = J2{-iub) n F^' n \r). (3.17) 



n=0 
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Then we can indeed write G^ which gives the entire first order correction to the metric 



compactly in Fourier space as 



<%\r, U ) = #F^r,u>)n ij (u>). (3.18) 



This equation is analogous to ( p. 23 ) where the full correction to the metric at the first 



order has been captured by a single function f^\r,u) after Fourier transforming the time 
dependence in Fefferman-Graham coordinates. 

Finally, the equation for translation (|3.12|) when Fourier transformed in v, gives the 



following algebraic relationship between and 

F«(r, W ) = —^—J {1) (P = * (0) (rM ^ fe(0)( °- (3-19) 

If we were able to solve the ODE for / (1) (p, uj) given by ( |2.33| ) exactly, we could have easily 
used the algebraic equation above to know F^(r,u) exactly. However, this has not been 
possible, so we have taken care of the dependence on u, or equivalently summed the time 



derivatives recursively. We have already translated the same strategy here in (3.14j) 



3.3 Translation at the second order in the amplitude expansion 

As in the case of the first order in the amplitude expansion, it is instructive to begin with 
how the changes in coordinates ( |3.8|) take their form at the second order in the amplitude 
expansion. Taking into account manifest invariance under spatial translations and rotations, 
clearly <3>( 2 ) (r, v) and (r, v) specifying the change of coordinates at the second order in 
the amplitude expansion, should take the form, 

oo n n , , s a / j \ b 

$ (2) M) = E E & 9+ "$ (2 '"' m) (r) E (j-j ^ v At) 

n=0 m=0 a,b=0 ^ ^ ^ ' 

n+m is even a+b=n 

\a—b\=m 

oo n n , , v a / , \ b 

k W(r,v) = E E b^k^(r) E £ 7r " (v) ' (3 - 20) 

n=0 m=0 a,b=0 ^ ' ^ ' 

n+m is even a +b=n 

\a—b\=m 

where, by definition, $( 2 ' n > m )(r) and k^ 2,n,m \r) are dimensionless functions, and n + m is 
even ensures that the sum over m is over even integers if n is even, and is over odd integers 
if n is odd, so that both (n + m)/2 and (n — m)/2, i.e. a and b, are integers. 

In Fourier space, we can put these changes of coordinates, compactly in terms of just two 
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functions, $( 2 )(r, cj, oji) and k^ 2 \r, oj } wi), through the following relation, 

/DO /'OG 
/ dwj $ {2) (r, w,Wi)7r M (wi)7r M (a; - u^), 
-OO «/ — OG 

/*oo /*oo 

A;( 2 )( r ,t;) = b 9 / du e~ iuJV dui k {2 \r, u, wi)7r M (wi)7r M (w - Wi), (3.21) 



where, 

^(r.w.wO = E E (-*) n K T ( W - W i)^ 



oo n 



n=0 m=0 

n+m is even 



n — m n4-m / \ 

+w~(w - wi)^)$ (2 ' n ' m) (r), 



oo n 



n=0 m=0 

n.+m is enen 



^ )^ m )(r). (3.22) 

We note that, by definition $^(r, w, Wi) and fc^(r, a>,a>i) are dimensionless functions. 

We can readily obtain the equations for & 2 \r, u, ui) by expanding the equations for 
change of coordinates (|3.9| ) up to second order in the amplitude expansion, and then Fourier 
transforming the time dependence. Thus we have, 

d$(0)( r ) 0$(2) dfc (0)( r ) ^(2) 



+ - 



1 / ^ ( (0) 



00 



2 \ dp I \ dr 



(0) _ *fol^N dk(2) f- . \ ,,,.rf$ (0) (r) 



»S , (p = * (0) (r))^-(r )W)Wl ) (iw 



<9r c?r 



(p = $(°)(r)) 



dp J dr 
+ 2$(°)( r )) $( 2 )(r,w,a;i) 

-iurf ( P = $(°)(r))^lA;( 2 )(r, W ,a; 1 ) 

= -i/f (P = *<°>(r),",«i) e-^ (0) «^ . (3.23) 
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The solutions are uniquely specified by the boundary conditions, 



$ (2) (r, u, ui) 



O 

o 



b 9 r 9 atr^oo, 



at r — > oo. 



b 9 r 9 



(3.24) 



These boundary conditions follows from simple dimensional analysis. The leading behavior 
of $(°)(r) and k^(r) at the boundary r = oo is 1/r, and is independent of b. The leading 
bahavior of $( 2 )(r, v) and k^(r,v) also should be independent of b, so dimensional analysis 
shows that their leading behavior should be 1/r 9 . This begets the above boundary conditions. 

The eqs. ( |3.23| ) which are ordinary differential equations in r are unfortunately not 
solvable for arbitrary u and U\. Therefore, we can adopt the same strategy as in the previous 
section, to solve these exactly first when u = uj\ = 0, so that we know $( 2 > '°)(r) and k^ 2 '°'°\r) 
exactly and then sum over the dependence over u> and U\ in ( |3.22| ) by obtaining $( 2 ' n - m )(r) 
and fc( 2 ' n ' m )(r) recursively. We will not repeat this here because it will turn out that these 
exact recursion relations will not be important for the regularity analysis. 

We also note that the eqs. ( |3.23| ) are first order ordinary differential equations in r. 
Given that the equations for translation of the metric which we will derive next are algebraic, 
we bypass solving second order differential equations with constraints, involved in directly 
obtaining the metric in Eddington-Finkelstein coordinates from Einstein's equation. 

The G^v in the ingoing Eddington-Finkelstein metric ( [3.6|) , similarly, should take the 
following form at the second order in the amplitude expansion, 



G$!(r,v) 



n=0 m=0 

n+m is even 



6+n Ei(2,n,m) . 



E 

a,b=0 
a+b=n 
\a— b\=m 



n=0 m=0 

n+m is even 



6+n 77.(2, n,m) . 



^ E {£) «m (i ) 



+E E b6+np i 



6+n 77>(2,n,m) . 



n=0 m=0 

n+m is even 



a,b=0 
a+b=n 
\a— b\=m 



E 

a,b=0 
a+b=n 
\a— b\=m 



d y 



-Sij ( — 1 n,„,{ '')( — ) 7Y lia ( r 



dv J 



dv ) 



. (3.25) 



By definition, F^ 2 ' n ' m ^s are dimensionless functions. Once again the sum over m runs over 
even intergers if n is even, and over odd integers if n is odd. 
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(2) 

After Fourier transforming the time dependence, in analogy with g^J (p, 00, oo\) in the 
Fefferman- Graham metric at the second order in the amplitude expansion given by ( |2.27|) , 

(2) (2) 

G}J(r, co, uji) can be compactly written in terms of three dimensionless functions (r, oo, ooi), 
for % = 1, 2, 3 as 



r*oo />oo 

G oo( r > v ) = b 6 du e~ luJV I dio 1 F^ (r,uo,oo 1 )-K pq {oo 1 )'K pq {oo - u)i), 



-00 j —00 

OO /"OO 



GfMr,v) = 6 6 / doo e'™ I duo, F?\r,u, 00^(00^(00 - Ul ) 



00 j —00 

OO POD 



+ y / due l0JV I d0Oi F^\r,0O,0Oi)[lT ik (0Oi)TT kj (0O - 00!) 
J —00 J —00 

+7l ik (00 - Wl)7Tfcj(u;i) 

2 

—^8 ij 'K pq (uo 1 )'Kp q {uo -oox)], (3.26) 

where 

00 n 

^ 2) (r )W)Wl ) -EE (^"(^^-^^ 

n=0 m=0 

n+m is enen 

+w 1 ^(a;-a;i) Ii ^)F i (2 ' n ' m) (r),for 2 = 1,2,3. (3.27) 

The equations for the translation of the metric (ft.lOD at the second order in the amplitude 
expansion, after Fourier transforming the time dependence implies, 

if } (Wi) = — |^(/s (2) (P = $ (0) (r),c,ci)e-^ (0) « 

+b (^-j (p = ^(r))^ 2 \r,oo,oo 1 ) 

—i2bug^k^ 2 \r, 00, 001) 
-26Gg ) (0^ (0) W* (2) (^^^i)) ) 
ifV,"^) = — |— .(/f (p = ^(r),c,ci)e- fc(0) « 



6 3 

-2for 2 $ (Q) (r)$ (2) (r,a;,a; 1 )), 

F< 2 >(r, W>Wl ) = — *fV( p = *V\ r ),u,, Ul )e^<r\ (3.28) 
(q>W(rJJ 

(2), 



Here we have used the explicit second order form of g^J(p,u,ooi) as given by fl2.27|) . The 
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above equations are algebraic. We recall that ( |2.41| ), ( |2.46| ) and fl2.49| ) gives the equations of 

(2) 

motion for fl (p, oj, ui), for i — 1, 2 and 3 respectively, and they have unique solutions given 
the boundary conditions ( |2.52|) and (|2.53|) . If we use these and the solution for $( 2 )(r, u, Ui), 



we obtain the full translation of the metric into the ingoing Eddington-Finkelstein coordi- 
nates through these algebraic equations. 

In practice, we need to solve for the time dependence, or equivalently the u and uj\ 
dependence of F^ 2 \r, oj, Ui) for % = 1,2,3 given by ( |3.27| ), by solving the eqs. ( |3.2£| ) for 
to = u\ = first, and then obtaining F^ 2,n ' m \r) and k^ 2 ' n,m \r) recursively. We need to use 
the recursion relations for y.( 2 ' 2ri ' 2m ) (p) given by (|2.57| ), (|2.67|) and (|2.68| ) for i = 1,2 and 3 
respectively and also the recursion relations for $( 2 ' n > m )(r). 

It will turn out that for the regularity analysis, we will only need recursion relations for 
(r) which, according to (|3.9| ), can be obtained without knowing $( 2 ' n > m )(r). We get 
the recursion relations, 

t.2 2 / m 

if ■ B ' m) (r, u , Ul ) = ff Y y f^ m ^ (p = $(°) ( r )) 



+ j2 f { i' 2p,m+q) (p = ^ (0) w) ] 1 j 

9=0 

when n is even, 



b n " 2 Pq\(n -2p-q)\ 



n-l 
U2 2 



2p—m \ h (0) { \\ n ~ 2p 

+ J2f?> 2M) (p = &°\r)) 

g=0 



(jfe(0)( r ))' 



b n - 2 Pq\(n -2p-q)\ 
when n is odd, (3.29) 



where in turn f[ 2 ' r ' s \p) are defined through the recursion relations ( [2.57 ), and vanish unless 



both r and s are even by their definitions in ( f2.28| ). As a specific instance, 

F p.o,») (r) = *V ^ og ^_±_^\ (3 . 30) 

(2) 

which is the only contribution to F± if 7T»j is time- independent. 

Clearly, we can use our method for translating the metric into ingoing Eddington- 
Finkelstein coordinates from the Fefferman-Graham coordinates at higher orders in the 
amplitude expansion, summing over all time derivatives at each order. 
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4 Regularity analysis of the metric and new phenomeno- 
logical parameters 

As discussed in the Introduction, it can be expected that the metric will be manifestly 
regular at the future horizon in an appropriate coordinate system order by order in the 
amplitude expansion provided we sum over all time derivatives. The amplitude of the non- 
hydrodynamic shear-stress tensor becomes arbitrarily small when compared to the pressure 
close to equilibrium. However, its time derivatives are not small compared to the temperature 
even close to equilibrium. This necessitates summing over all time derivatives while treating 
the amplitude expansion perturbatively. 

When the energy-momentum tensor is purely hydrodynamic, the derivatives of the hydro- 
dynamic variables are small compared to the local temperature. The derivative expansion, 
therefore, can be employed perturbatively. This logic was instrumental in the construction 
of the metric perturbatively in the derivative expansion, which was regular at the future 
horizon for appropriate choice of the transport coefficients. We expect that in our case, sim- 
ilarly, we will be able to extract an equation of motion for the shear-stress tensor order by 
order in the amplitude expansion, so that when the shear-stress tensor follows this equation 
of motion with appropriate values of the phenomenological parameters, the metric will be 
regular order by order in the amplitude expansion. In this Section, we will outline the details 
of this analysis. 

Before we go through the details, we need to consider if we can do the regularity analysis 
more efficiently by calculating curvature invariants like R fiupcr R^pa at the horizon. Unlike 
the hydrodynamic case, where the regularity simply involves the transport coefficients, and 
regularity analysis simply tells us we have a regular horizon provided some algebraic equa- 
tions involving the transport coefficients are satisfied, we expect that the regularity condition 
here will involve an equation of motion for TTij(t). Calculation of curvature invariants, which 
we will not reproduce here, also indicates so. However, a curvature invariant can at best 
reproduce the trace of the equation of motion for iiij(t), with for instance, a linear combi- 
nation of vrjj(t) and its time-derivatives. Therefore, one needs to compare several curvature 
invariants, to extract the equation of motion itself. This has proved to be a very difficult 
task when we have tried to do it. More than this practical issue, one cannot be sure of 
the regularity unless one has compared all possible curvature invariants and this requires 
utilising geometric properties of the spacetime. It seems that it is very difficult to do this 
also. The most elegant way to do the regularity analysis, it appears, is to go to a coordinate 
system where the metric can be argued to be manifestly regular. Here, we will show that we 
can indeed do so in the ingoing Eddington-Finkelstein coordinates. 

We will first do the regularity analysis at the linear order, that is at the first order in the 
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amplitude expansion. Then we will extend our analysis to higher orders in the amplitude 
expansion. This will allow us to extract the equation of motion for TTij(t) in the amplitude 
expansion, such that when this is obeyed, the corresponding metric in the bulk has a regular 
future horizon. Finally, we will connect our results with our previous work and see how we 
can fix some more general phenomenological coefficients at strong coupling. 



4.1 Regularity at the linear order 

At the level of linearized gravity perturbations with right asymptotic behavior, quasinormal 
modes are those which are ingoing at the future horizon. Though there is no rigorous proof, 



there are strong arguments p6| [30| to suggest that for linearized gravity perturbations, the 
ingoing boundary condition at the future horizon is equivalent to manifest regularity of the 
metric at the future horizon in ingoing Eddington-Finkelstein coordinates. In other words, 
when the metric is perturbed with the quasinormal modes, not only combinations of the 
linear perturbations of the metric invariant under linearized diffeomorphisms, but the full 
metric tensor itself will still be manifestly regular at the future horizon in ingoing Eddington- 
Finkelstein coordinates. 

At the first order in amplitude expansion, the dependence on the perturbation of the non- 
hydrodynamic shear-stress tensor is linear, so the manifest regularity at the future horizon 
in ingoing Eddington-Finkelstein coordinates should require that the shear-stress tensor has 
time dependence given by a linear superposition of the quasinormal mode frequencies at zero 
wavelengths. 

In the computations of the quasinormal modes for AdS black branes, usually the ingoing 
boundary condition is imposed first and the dispersion relation is then extracted by requiring 
that the Dirichlet boundary condition is obeyed at asymptotic spatial infinity. This would 
imply that the non-normalizable mode of the perturbation vanishes. However, in the ap- 
proach employed here, we first find the metric corresponding to an arbitrary 7Ty(t) in a flat 
background in the dual theory, thus fixing both the non-normalizable and the normalizable 
modes of the metric perturbations. These determine the solution uniquely. The quasinormal 
modes can be expected to be extracted by requiring manifest regularity of the metric at 
the future horizon in ingoing Eddington-Finkelstein coordinates. Our method can be read- 
ily generalized at the non-linear level, besides it will give us the most general metrics with 
regular future horizons dual to homogeneous relaxation. 

To see how this works, it will be instructive to review the case of the purely hydrodynamic 
shear-stress tensor in this new light. 

Let us consider a fluid at constant four- velocity w^ ^ and constant temperature T^°\ It 
is convenient to define as as before, so thatl/M ) is the radial location of the 

future horizon in the dual geometry in ingoing Eddington-Finkelstein coordinates. Further 
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we define the projection tensor in the spatial plane orthogonal to vf®" as 



P (oW = u (o)n u (o)u + v nu_ ^ 

The most general purely hydro dynamic shear-stress tensor linear in the fluctuation in veloc- 
ity, 5u^ and the temperature, 5T, is 



26( ) 3 

a ' 



46(o) 4 Mop [j/oj J 

<x^ + 0(e 2 ), with (4.2) 



where 0(e 2 ) denote terms at second and higher orders in the derivative expansion. Also 7 
is a dimensionless quantity defined through 

7 = ^, (4.4) 
s 

with i] being the shear viscosity and s the entropy density. 

The corresponding metric in ingoing Eddington-Finkelstein coordinates up to first order 
in the derivative expansion is, 

ds 2 = -2{uf + 5u^)dx»dr + G^cbfdx", (4.5) 

G = r 2 P (0) + (-r 2 + 1 W (( V 0) — ( —\ M (<V 0) 



^ [iv 1 V ' ^(0)4^2 ' A 4 ^ 6(0)4 r 2 I £,(0) i A 4 ^ 

-r ^(u (0)Q 9 Q )^ + uf(u^ a d a )5u^ - ^ufu^\d a 5u c 
+ ( 26 < Vlf ( 6 W>r) + ^6<V log (l - ^ ) <v + O (e 2 ) , 

where defined in ( f4.2| ), 0(e 2 ) denotes terms at higher orders in the derivative 

expansion, and 



H(x) = i (log f ( x+1 )y 2 + 1 ) ) _ 2 M:cUlu[x] - ) . ( Ui) 



The metric above is just the linearization of the metric for arbitrary 7 in ]I4|], which gener- 
alizes the metric in [[L0] for 7 = 1. 
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We see that corrections to G^ in the metric in ingoing Eddington-Finkelstein coordinates, 
which is linear in 5u^ and 5T, can be decomposed into four categories based on whether the 
components are in the direction of the flow 

u (o)p or 

orthogonal to it. The first category is 
longitudinal-longitudinal, i.e. proportional to the second category is longitudinal- 

transverse , i.e proportional to terms like u^5u u , or u^\u^ ■ d)5u u ; the third category is 
pure trace and transverse-transverse, i.e. proportional to Pjh); and the fourth one is traceless 
and transverse-transverse, i.e proportional to terms like a^. This division can obviously be 
done at any order in the derivative expansion. 

We also see that, at the first order in the derivative expansion, the leading divergence 
at the late-time horizon which coincides with the unperturbed horizon r = is ((7 — 

l)/4) log(&(°V — 1)0^ and belongs to the fourth category, which is traceless and transverse- 
transverse. This divergence vanishes provided 7 = 1, i.e. when r]/s — 1/Att. The other 
divergences are of the same category, and are coefficients of — 1) log(&(°V — 1) and 

(p(°)r — l) 2 log(&(°V — 1), both of which are proportional to (7 — 1)<t^, which vanish for the 
same choice of rj/s. These terms are potentially divergent at the horizon because their first 
or second derivatives, or both, with respect to r diverges at the horizon, i.e. at r = l/b^°\ 

If we do not fix 7 (or equivalently rj/s) at the first order and obtain the metric at higher 
orders in the derivative expansion, we find various other kinds of divergence in all four 
categories of G^ v in the ingoing Eddington-Finkelstein metric. When 7 = 1 all divergences 
in all four categories of G^ v , with higher powers of the logarithm, like (log(£/°V — l)) 2 and 
proportional to derivatives of a^ v vanish. However, only three divergences of lower orders 
remain. 

These are as follows, one is a log(6^°V — 1) divergence, another is a (&(°V — 1) log(&(°V— 1) 
divergence and the other one is a — l) 2 log(6^°V — 1), all of which at the second order 
in the derivative expansion, for instance, are proportional to (u^ ■ d)a lxv (which is Weyl 
covariant at linear order in the velocity fluctuation Su), and belongs to the traceless and 
transverse-transverse category as it is orthogonal to u^. The cancellation of these divergences 
requires fixing the only linear transport coefficient at the second order in the derivative 
expansion in a flat background, whose contribution to the shear-stress tensor at the second 
order in the derivative expansion, is proportional to (u^ ■d)a flI/ . Further, all these divergences 
are proportional to a polynomial which is linear in the arbitrary choice of this second order 
transport coefficient and vanish when this transport coefficient takes the same value, which 
is (2-log2)/6(°) 2 . 

It can be shown that the feature just mentioned persists up to arbitrary higher orders 
in the derivative expansion. The transport coefficients at any order can always be fixed by 
taking care of the log(6^°V — 1) divergence in the traceless and transverse-transverse part 
of G^ u , and once these are fixed at any given order, there is no other kind of divergence in 
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any category of G^ v provided the transport coefficients at all lower orders have been chosen 
correctly. Also this log(ir°V — 1) divergence is linear in the arbitrary choices of the transport 
coefficients at that order. To summarize, the cancellation of the log(6^°V — 1) divergence 
ensures that all other divergences vanish order by order in the derivative expansion, but not 
vice- versa. 

Therefore, the regularity condition at the future horizon is simply that the coefficient 
of log(fe(°V — 1) term in the traceless and transverse-transverse part of should vanish, 
and this fixes the value of all transport coefficients to all orders in the derivative expansion. 
Up to second order in the derivative expansion, explicit calculations have shown that the 
values of transport coefficients determined by regularity at the future horizon in ingoing 
Eddington-Finkelstein coordinates match with those obtained from the dispersion relations 
of the hydrodynamic shear and sound quasinormal modes. 

Even in the case of homogeneous relaxation, the regularity condition at the horizon should 
be the same. The argument for it is easy. The regularity condition should be a condition 
on 7Tij{t), therefore it is independent of r. If the divergence at the late-time horizon has 
to vanish at all times, the coefficient of each singular term in the series expansion of G^ v 
in r has to vanish at r = 1/6, where the future horizon is located at late time. f\ The 
singular terms should also appear in the traceless and transverse-transverse part of G^, 
since all the corrections to the metric at the first-order in the amplitude expansion being 
porportional to 7Ty(t) and its time derivatives as shown in (|3.12| ) and (|3.13|) , are traceless 
and transverse-transverse. 

Further, the singular term in the series expansion of G^ v in r at r = 1/6, whose coefficient 
should give a sufficient condition for all other divergent terms in r to vanish at the late time 
horizon, can be identified as the log(6r — 1) term, by looking at the case when TTij is time- 
independent. In this case, the entire contribution to the metric at the first order in the 
amplitude expansion, according to ( |3.13| ) and (|3.15|) is given by 



Gg>( r ) = b 2 F^\r)^ = -6V log (l - -L) *■<,. (4.7) 

Indeed the vanishing of the log(6r — 1) term, which is the leading divergence at r = 1/6 here, 
gives the sufficient condition tt^- = 0, for regularity. This is expected because the black hole 
cannot support any tensor hair. 

Time- dependence of ^(i) brings in divergences with arbitrary higher powers of the 
logarithm, i.e. terms like (log(6r — l)) n , where n > 1 is an integer; however it is the time- 
dependent corrections to the coefficient of log(6r — 1) that should control all the other 
divergences. In fact, simple inspection of FW(r, to) through (|3.19|) , shows that the coefficient 



3 We recall that b is constant here to all orders in the amplitude expansion. 
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of log(6r — l) n term for n > 2 is 0(u n ). As these coefficients of such singular terms have 
higher powers of u at the leading order itself, these should be total time derivatives of the 
coefficient of the log(6r — 1) term, which is 0(1) in the u expansion. Thus, the vanishing of 
the latter can ensure that other divergences vanish, but not vice-versa. 

Therefore, the regularity condition at the future horizon is that the coefficient of the 
log(6r— 1) term in the traceless and transverse-transverse part ofG^ in the ingoing Eddington- 
Finkelstein metric should vanish, which implies that ^(t) should obey a linear differential 
equation that can be calculated from Einstein's equation. 

If the perturbations corresponding to quasinormal modes are indeed manifestly regular in 
the ingoing Eddington-Finkelstein coordinates, by the argument above, the only possibility 
for the regularity condition at the horizon is that the coefficient of the log(6r — 1) term should 
vanish, so that it will also provide a sufficient condition for all other divergences to vanish. 

We need to understand now how we can reproduce the quasinormal mode frequencies at 
zero wavelengths from this regularity condition. Let the coefficient of the log(6r — 1) term 
in the series expansion in r of F^(r,u>) about r = 1/6 be D^\u). Then our regularity 
condition at the linear order is 



where the integration is done with the contour prescription mentioned before, which picks up 
contributions only from the poles of the analytic continuation of D^\uj) 11.^(00) as a function 
of the complex variable oj in the lower half plane. By definition, D^\co) is a dimensionless 
function. We have also observed before that TTij(u)) is identical if it is the Fourier transform 
of iTij(t) or 7Tij(v), because the latter are the same functions, though of different variables. 
So, the regularity condition on 7iij(uj) leads directly to the regularity condition on 7Ty(t) after 
doing Fourier transform in t. 

Further, near uj = 0, we can write, 



where clearly is the coefficient of log(6r — 1) term in the series expansion of F^' n \r) 




(4.8) 




(4.9) 



n=0 



about r = 1/6 as is obvious from ( 3.17|) . By definition, D^ 1 '^ 7 



s are dimensionless numbers 



which can be explicitly evaluated using (|3.14j) . For instance, 



D 



"I, 




etc. 



(4.10) 



However, this power series ( |4.9j ) can have a radius of convergence. So, it is necessary to 
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define an analytic continuation of D R '(to) beyond this radius of convergence in the lower half 
plane, for doing the Fourier transform using the mentioned contour prescription. It is also 
important to realize that it is the analytic continuation of (oj) and 7Ty(u;) in the lower 
half plane as a function of the complex variable u which is required to make connection with 
the quasinormal modes. 

If the analytic continuation of D^\oo) have simple zeroes at ± Lo^ n —iujj n , with ujj n > 0, 
in the lower half plane, and ffij(oj) have at most simple poles at ± u>R n —iujj n in the lower half 
plane, then it follows that D^\u)TCij(uj) is analytic in the lower-half plane and its Fourier 
transform vanishes by our contour prescription as required by the regularity condition (|4.8|) . 
Thus the regularity condition implies that in the lower half complex u plane, ^ (cu) should 
have the following general form 

M») - 1 {I) ( Mt .-,V-^ ) + „,-l\ u J + ^ 

where f(u) for large values of \cu\ decays like 0(l/\ui\ 2 ). So for t > 0, the Fourier transform 
of 7Tjj (co) is 

oo 
n=0 

which will match with the prediction of the quasinormal modes if ± ujR n — itoi n , with ui n > 
are the quasinormal mode frequencies in the scalar channel when the wavelengths vanish. 
This spectrum in the scalar channel has been numerically calculated in 

The coefficients of other singular terms, like log(6r — l) 2 , for instance, in in ingoing 
Eddington-Finkelstein metric at this order, should also take the form, 

/oo 
duj e- iu}V D^(cu)7r i:j (uj) = 0, (4.13) 
-oo 

where D^(co) is the coefficient of the log(6r — l) 2 term in the of F^(r,u) obtained in the 
previous section. We require that the analytic continuation of D^(uj) also should have 
simple zeroes at the same points ± u>R n — iuj n in the lower half complex plane. We recall 
that, in the purely hydrodynamic case, for instance, at the second order in the derivative 
expansion, the analogous coefficient was a second order polynomial in 7 (i.e. Atttj/s), which 
had a simple zero at 7 = 1, just like the coefficient of the log(6r — 1) term. Thus the regularity 
condition (|4.8| ) should provide sufficient condition for cancellation of all divergent terms. 

To summarize, we see that at the first order in the amplitude expansion, regularity 
analysis at the late-time horizon simply implies that the metric perturbations are composed 

10 In it was also found that the quasinomal spectrum at zero wavelengths in the scalar channel is well 
approximated by 2nTn (±1 — i) for large n. 
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of the infinite tower of quasinormal modes in the scalar channel with zero wavelengths. 

To check this explicitly, we need to find the power series ( USD defining D%\u) about 
u = explicitly. We have been able to define this implicitly as are the coefficients of 

the \og(br — 1) terms of F^' n \r), which have been obtained through the recursion relation 
( P-14Q . If we can find an explicit recursion relation in terms of D^' n ^ directly, we will be able 



to see how we can analytically continue the power series (4.9) in the lower half complex oj 
plane, and check explicitly if it has only simple zeroes and if they match with the quasinormal 
spectrum. This explicit recursion relation is unfortunately a hard combinatorial problem and 
we have been unable to solve it here. 

4.2 The full regularity analysis 

When the shear-stress tensor is purely hydrodynamic, the metric continues to be manifestly 
regular at the late-time horizon in the ingoing Eddington-Finkelstein coordinates at the non- 
linear level too, order by order in the derivative expansion. It is this feature that allows us to 
determine the non-linear transport coefficients which do not affect the dispersion relations of 
the hydrodynamic shear mode and sound modes, order by order in the derivative expansion. 
At second order in the derivative expansion, in fact, we have three such non-linear transport 
coefficients, whose values have been determined by this method [[Hj . 

An arbitrary quasinormal mode, even if it is non-hydrodynamic, can be expected to 
be manifestly regular at the future horizon in ingoing Eddington-Finkelstein coordinates. 
However, at the non- linear level, when we build solutions corresponding to a "liquid" of such 
fluctuations, we have no strong general argument why the solution should be manifestly 
regular at the late-time horizon in ingoing Eddington-Finkelstein coordinates. 

Fortunately, in our present case of homogeneous relaxation we can strongly argue that 
the regularity of the metric at the late-time horizon will still be manifest in the ingoing 
Eddington-Finkelstein coordinates, order by order in the amplitude expansion. 

The metric in the Fefferman- Graham coordinates always maintains the symmetries of 
the dual boundary configuration, in this case the symmetry under spatial translations and 
rotations. However, this metric is not manifestly regular at the late-time horizon even 
when the regularity condition is satisfied. The change of coordinates which will achieve the 
manifest regularity at the late-time horizon can also be expected to respect the symmetries of 
the dual boundary configuration. This means that the spatial Fefferman-Graham coordinates 
z*'s should translate to the new spatial coordinates x l 's in order to maintain manifest spatial 
translation and rotational invariance up to rotations and shifts which depend on the new 
radial coordinate r and new time coordinate v only. So, 

z i = ifyr, v)x j + S*(r, v). (4.14) 
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However, it is not possible to write an orthogonal matrix R l - and a shift vector S l using just 
iTij and its derivatives. The rotational matrix and the shift vector, therefore are constants, 
in which case of course the metric remains invariant. 

Further, to maintain manifest spatial translational and rotational symmetries, we obvi- 
ously require that the Fefferman-Graham radial coordinate p and time-coordinate t, should 
be translated to new radial coordinate r and time coordinate v without involving the spatial 
coodinates, so that both p and t will be functions of r and v only. Besides, we want manifest 
regularity only at the future horizon and not at the past horizon. So the obvious choice is the 
ingoing Eddington-Finkelstein gauge where g rv — g vr — 1 (rather than -1) and g rr = g ri = 0, 
as at the zeroth order this choice makes the regularity at the future horizon manifest. We 
can expect this choice to achieve manifest regularity at the future horizon to all orders in 
the amplitude expansion. 

Repeating the analysis at the linear level or at the first order in the amplitude expansion, 
the regularity should be guaranteed by the vanishing of a certain singular term in the series 
expansion in r of G^ in the ingoing Eddington-Finkelstein coordinates, otherwise we cannot 
guarantee regularity at all times. Further, the choice of this singular term can be made by 
doing the regularity analysis at the linear level, so the regularity condition up to any given 
order in the amplitude expansion is the vanishing of the coefficient of the log(6r — 1) term 
in the traceless and transverse-transverse part of in the ingoing Eddington-Finkelstein 
metric. 

Clearly, the regularity condition up to second order in the amplitude expansion can be 
written in the form, 

/oo 
du e-™ Z^H^M 
-oo 

6 6 f°° r 
+— / duiD^\u,Ui) (Ti ik (ui)'n kj {uj-uj l ) + Ti ik (uj-uj l )'n kj {uJi) 

J -oo L 



-g5tf7r„(wi)7r„(w-wi))] + 0(5 3 )] =0, (4.15) 



where all the integrations are done with our contour prescription which picks up contributions 
only from poles in the lower half plane. Further, D^\u,uji) is dimensionless like D^\co) 
and can be readily identified as the coefficient of the log(6r — 1) term in the series expansion 

(2) 

of F{ '(r,uj,uji) at the late-time horizon r = 1/6. 
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We can also series expand D R '(u), U\) near uj = Ui = as below, 



r>£V,«i) = E E (-muTiu-ui)*? 



n=0 m=0 

n+va is even 



(4.16) 



where .D 



(2,n,m) 



are dimensionless numbers and the sum over m runs over even integers only if 



n is even, and odd integers only if n is odd, so that both (n+m)/2 and (n—m) /2 are integers. 



Comparing with (|3.27|) and recalling that D R (lu, u\) is the coefficient of the log(6r — 1) term 
of the series expansion of F^ 2 \r, oo, oji) at r = 1/6, we obtain that D R ' n ' m ^ is the coefficient 
of the log(6r — 1) term of the series expansion of F^ 2 ' n ' m '(r) at r = 1/6. For instance, 



D 



(2,0,0) 
R 



1/2, etc. 



(4.17) 



Given the explicit recursion relations for F{ 



(2,n,m) , 



)s given by ( |3.29| ) and ( |2.57| ), we in- 



(2) 

deed obtain recursion relations for D R (u,ui) implicitly. A more explicit answer, which is 
essentially a combinatorial question, will tell us how we can analytically continue the series 
( |4.16| ) in the lower half complex u and u% planes, to obtain the regularity condition explicitly 
through the Fourier transform ( |4.15| ). Unfortunately, we have not been able to achieve this 
here. 

However, the series expansion ( |4.16|) allows us to write the regularity condition fl4.15|) up 
to second order in the amplitude expansion in the following form, 



Csij (t) = ^E^g)" 

OO TL Th 

+^ 6 E E Di n n ' m)bn E 



n=0 m=0 

n+m is even 



a,b=0 
a+b=m 
\a—b\=m 



It nk \dt> ^ 



hJ dY 

3 



5 ^Jt) npq 



+ O(S 3 )=0, {A.li 



where again the sum over m is over even integers only when n is even, and over odd integers 
only if n is odd, so that both a and b are integers. Further, the reader may recall that the 
above equation has been derived in the global comoving frame where the velocity of the flow 
vanishes. However, this equation can be also written in a Lorentz-covariant form, with t 



replaced by —u.x and d/dt replaced by u ■ d. Also the covariant form of Hi 



IS 7T 



, which by 



definition is such that its non-vanishing part is the projection to the spatial plane orthogonal 
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to u p , i.e. P fl p P u a ir pr7 , with P^ p being the spatial hyperplane projection tensor u p u p + 5 p p . 
The regularity condition in gravity stated here, which is the vanishing of the log(6r — 1) term 
in the series expansion of traceless and transverse part of G^, is also Lorentz-covariant and 
therefore automatically yields this form. 

The non-linear terms in this equation of motion for 7Ty(t) can have drastic effects. A 
specific question can be elegantly framed in the Fourier space. We know from the discussion 
in the previous subsection that the regularity condition at the linear order implies that 
71^(00), when analytically continued as a function of u in the LHP, should have an infinite 
series of discrete poles implying exponential decay in real time. The question is how the 
behavior is modified at the non-linear level, do the residues of the poles get related and the 
subleading 0(l/u) 2 ) asymptotic behavior constrained by the non-linear corrections, or is the 
modification more radical, like appearance of branch cuts, etc.? 

We conclude our regularity analysis here by emphasizing again that the eq ( 4.18[ ) is valid 
only when the amplitude of the non-hydrodynamic is small compared to the pressure, 
though its time derivatives can be large. Further, such solutions should be smoothly con- 
nected to the equilibrium solution, where 7Ty = 0, because the dual metric in gravity from 
which this eq ( |4. 1 8| ) has been obtained, is connected to the equilibrium solution by smooth 
variation of the amplitude perturbation parameter. This also implies that the meaningful 
solutions equilibrate at late-time. It thus makes no sense, for instance, to look for time- 
independent solutions to ( |4.18| ) truncated up to second order in the amplitude expansion, 



though such non-trivial solutions discrete up to constant spatial rotations exist. Besides, one 
can also check that, for this instance, the amplitude is of 0(1) and hence the perturbation 
expansion can not be trusted anyway. 

4.3 Extracting new phenomenological parameters 

In the case of the Boltzmann equation, as discussed in the Introduction, we can construct 



conservative solutions for homogeneous relaxation, where we can obtain an equation of 
motion for the shear-stress tensor r K%j{t\ such that any solution of this equation lifts to a 
unique solution of the full Boltzmann equation. 

The shear-stress tensor is actually the local symmetric and traceless rank two veloc- 
ity moment tensor of f(x,v), the quasiparticle distribution function. In the homogeneous 
solutions of the Boltzmann equation, the hydrodynamic variables are constant in space and 
time. Further, in the conservative homogeneous solutions 7Ty(t) obeys an equation of motion 
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which can be obtained order by order in the amplitude expansion and takes the form 
dir. 



' 13 



dt 



BV\p, T)^ + 5 (4 ' 2 ' 2) (p, T) (n ik ir kj + 7i jk 7i ki - ^ fj 7r Iro 7r Iro ^ + 0(5 3 ), (4.19) 



where and _g( 4 > 2 ' 2 ) can be calculated from the collision kernel of the Boltzmann equation 
as functions of the thermodynamic variables, the density p and the temperature T, and 
0(5 3 ) denote terms third order in the amplitude expansion and beyond fll5 |. We recall that 



the velocity is constant, and in a Lorentz covariant notation, t = —u ■ x and d/dt — u ■ d. 
Also the covariant form of tt^- is the non- vanishing part of n^, which is the projection to 
the spatial plane orthogonal to u p , i.e. P^ p P u u tt P(T , with P^ p being the projection tensor 
u^u p + 5 p . In the conservative solutions, all velocity moments of f(x,v) other than the 
constant hydrodynamic variables are algebraic functions of 7^ and their time derivatives, so 
they do not have independent dynamical parts. 

The important difference of ( |4.19| ) derived via the Boltzmann equation which holds at 



small 't Hooft coupling, from ( 4.18 ) which holds at large 't Hooft coupling A, is this equation 
is only first-order in time. In fact, we know that such a Boltzmann equation is equivalent to 
perturbative non-Abelian gauge theory |TJ| |T7| at sufficiently high temperature, so this is in- 
deed a feature of perturbative non-equilibrium non-Abelian gauge theories. The higher order 



time derivatives in ( 4.18|) therefore have their origins in the intrinsic energy-time uncertainty 



of quantum dynamics going beyond the semiclassical approximation in the Boltzmann equa- 
tion, according to which we cannot define the energy of the excitation at short time scales. 
As of now, we do not know how to capture these effects in non-equilibrium non-Abelian 
gauge theories systematically. 

In case of a general inhomogeneous conservative solution, 7tij(x, t) can be split as ir^ h \x, t) + 
Tr^j\x,t). Here, 71^ is purely hydrodynamic and is given as an algebraic function of the hy- 
rodynamic variables and their spatial derivatives in the derivative expansion (in relativistic 
case, this is true only in a local inertial frame where the local mean velocity of the quasipar- 
ticles vanishes). On the other hand tt^J satisfies a equation of motion which is first order 
in time derivative (in the relativistic case only in the mentioned local inertial frame), and 
can be expanded in both the amplitude and derivative expansions. Both can be obtained 
systematically from the Boltzmann equation. In both cases, the derivative expansion refers 
to the spatial derivatives also in a local inertial frame where the local mean velocity vanishes. 

Because the Boltzmann equation itself is first-order in time, all higher time derivatives can 
be eliminated in favour of spatial derivatives in the conservative solutions, in an appropriate 



11 Here we assume that the Boltzmann equation is that of an underlying continuum system, it is not for 
instance a lattice Boltzmann equation. Further, the continuum system is not subject to external forces and 
its equilibrium state is isotropic. All these assumptions apply to gauge-theory plasmas. 
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local inertial frame just mentioned, so that the equation of motion of the hydro dynamic 
variables and the shear-stress tensor are first order in time. All other velocity moments 
of f(x,v) are algebraic functions of the hydrodynamic variables, the shear-stress tensor 
and their spatial derivatives, in the same local inertial frame, so do not have independent 
dynamical parts. 

The amplitude and derivative expansion become complicated if we go beyond the Boltz- 
mann approximation. The conservative solutions should also exist beyond this approxima- 
tion, as at large N and 't Hooft coupling A this would naturally explain why in the super- 
gravity approximation the universal sector of non-equilibrium states can be determined from 
the dynamics of energy-momentum tensor alone. This has already been discussed in the 
Introduction. However, there is no reason why the equations of motion for hydrodynamic 
variables and tt^ should be only first-order in time in a local inertial frame where the mean 
velocity vanishes. 

In presence of higher order time derivatives in the kinetic equation, we cannot treat 
the time derivatives of the non-hydrodynamic it^ perturbatively. Though tt}] 1 becomes 
arbitrarily small compared to the pressure close to equilibrium, its time derivatives do not 
become small with respect to average time between quasiparticle collisions even close to 
equilibrium, as it decays exponentially at linear order. This means that we should treat the 
derivatives of hydrodynamic variables perturbatively (including time derivatives), but we 
can treat the amplitude of 7rj"^ perturbatively only if we sum over all its time derivatives. 
However, close to equilibrium, we can consider a regime where the spatial derivatives are 
small. 

At large 't Hooft coupling A and large N, we can expect that the conservative solutions 
will form the universal sector of non-equilibrium states dual to the solutions of pure gravity 
with regular future horizons. In that case, we get the desired equation of motion of the 
energy-momentum tensor from the regularity condition on the late-time horizon in the bulk 
metric. In the purely hydrodynamic case before, and here in the case of homogeneous 
relaxation, such equations of motion have been thus obtained from gravity. However, we 
can also write down a phenomenological equation for the energy-momentum tensor which 
will hold for an abitrary configuration which equilibrates in the future. In absence of a 
known Boltzmann-like equation or any well-defined reliable formalism in the strongly coupled 
regime, the only guide to such a phenomenological equation is conformal covariance. 



In [15], we did propose such a phenomenological equation for the energy-momentum 
tensor, such that for the right values of the dimensionless phenomenological coefficients, it 
should reproduce all solutions of gravity with regular future horizons which equilibrate at 
late-time. However, we did not realize that we need to sum over all time derivatives of 
the non-hydrodynamic shear-stress tensor at each order in the amplitude expansion. This 
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did not allow us to compare with quasinormal modes correctly. We rectify this here and 
see how the analysis of homogeneous relaxation fixes some of the general phenomenological 
coefficients. 

To begin with, we will split the relativistic shear-stress tensor n^ u into a purely hy- 
drodynamic part tt^J algebraically determined by hydrodynamic variables and 7r^ whose 
evolution has dependence on both hydrodynamic and non-hydrodynamic parameters, and is 
therefore independent of the hydrodynamic variables dynamically. The amplitude expansion 
is concerned with only the non-hydrodynamic part, the dimensionless expansion parameter 
being its ratio with the pressure at late-time equilibrium. The hydrodynamic irffi is by 
definition an algebraic functional of the hydrodynamic variables and their derivatives (even 
time derivatives). The derivative expansion is concerned with spatial and temporal rate of 
variation of the hydrodynamic variables and only spatial rate of variations of tt^ , compared 
to the temperature at late-time equilibrium. Note the spatial variation of even 71"^™ is small 
close to equilibrium, but its time derivative is not. 

Time-derivatives involve (u ■ d), however if this acts on a Weyl-covariant quantity like 
\ the resulting quantity will not be Weyl-covariant. So, we need to define a Weyl- 



covariant form of the time derivative T>, following [37|, whose action on a Weyl-covariant 



tensor A^ turns out to be 

VA„ V = {u ■ d)A flv + -A^d ■ u) - (Afu v + Afu^) {u ■ d)u p . (4.20) 

Therefore, more precisely by time derivatives of 717^ we will denote successive actions of 
the Lorentz invariant and Weyl-covariant derivative T>. The spatial derivatives will denote 
actions of P u d v , where P u as before denotes u^u v + 8^ ', the projection on the spatial plane 
orthogonal to u^. 

The conservation of energy-momentum tensor implies 

V ( ^t^ + ttW + 7rg*>) = 0, (4.21) 

where b = 1/(ttT) as before, and all quantities are dependent on both space and time 
coordinates. We also use the earlier field definitions of T and such that (3/4) (vrT) 4 
is the local energy density and u^ 1 is the local four-velocity of energy-transport so that 
= 0. The traceless-ness of the energy- momentum tensor also implies that 
Tr(iv^) = Tr(ir( nh ^) = 0. With these conditions, 7r^ has five independent components 
only. As is an algebraic function of u M and T, and their derivatives, the nine independent 
variables parametrising the energy-momentum tensor are T, the three independent variables 
of the four velocity and the five independent components of 7T^ . 
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The form of the hydrodynamic part of the shear-stress tensor can be determined by Weyl 
covariance and the transport coefficients can be fixed by the regularity condition at the future 
horizon, order by order in the derivative expansion. Up to second order in the derivative 
expansion, we obtain 

(h) 2 2 -log 2^ 

2 ( 1 

+^r«"™ + v v a uap) + 0(e 3 ), (4.22) 

where 0(e 3 ) denotes terms higher order in the derivative expansion (where the expansion 
parameter e is the typical length scale of variation with respect to the temperature), and 

<V = \ p » ap J ( d aUp + dpu a ) - ^i%(<9 • u), 



1 



-P^P/(d a up-dpu a ). (4.23) 



The equation of motion of should be such that it is (a) Weyl-covariant and (b) we 
can consistently put n^J 1 ^ = 0, so that we can go to the pure hydrodynamic limit. In the 
case of the conservative solutions of the Boltzmann equation, this is always possible. Even 
in the case of gravity, we can construct solutions which are regular order by order in the 
derivative expansion when TijZ^ = 0, hence requirement (b) is essential. 

With these two requirements we can write the most general equation of motion for irfy!^ 



as 



E *™ = ^ + ^>a a , - \p^$**) (4-24) 



, A2 / 

2b ^ ' ^<wj 



IV V" n (2,n,m)J_ \ ^ 

n=0 m=0 a, 6=0 

n+m is even a+b=n 
\a—b\=m 



<r\a (nh)a<T\b (nh) 
jj, av 



+0(e 2 5, e5 2 ,5 3 ). 

Here, we note that the RHS is at least linear in the amplitude expansion (where the expansion 
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parameter 5 is the ratio of the typical non-hydrodynamic shear-stress with respect to the 
equilibrium pressure), so that Tr( nh > can be put to zero consistently. Further we have included 
terms only up to first order in derivative expansion on the RHS. The terms explicitly shown 
in RHS are thus 0(eS) and 0(5 2 ). We cannot include any term like <9 /i 7r^ h ' ) here because by 
( [4.21] ) it gets related to purely hydrodynamic quantities and their derivatives. By definition, 
all the phenomenological parameters Ai, A2, D^'^ and £)^ ,ra ' m ) are dimensionless. Further 
the condition n + m is even in the summation, denotes that sum over m is over odd integers 
only when n is odd and is over even integers only when n is even, so that both (n + m)/2 and 
(n — m)/2, i.e. a and b, are integers. This phenomenological equation is the right equation to 
consider so that we have indeed summed over all time derivatives of the non-hydrodynamic 
shear-stress tensor correctly. 

This equation ( |4.24|) , together with (|4.21|) and (|4.22|) give the complete equations of 
motion for the energy-momentum tensor in the general case. Further for the right choices of 
the phenomenological coefficients in ([4.24 ), the dual solution in gravity should have regular 
future horizons. 

We note that when and T are constants in space and time, and the fluid is at rest, 
the equation for conservation of energy and momentum given by ( |4.21|) implies that 7Tq^ = 
75 = ^io^ = 0' an d ^ij^ is a function of time only. Also, in such a case, the Weyl- 
covariant derivative T> reduces to an ordinary time derivative. We can readily see that the 
equation ( |4.24| ) rduces to ( [4.18| ), so that we recover the regularity condition at the late-time 
horizon for metrics dual to homogeneous relaxation. 

We cannot determine the phenomenological coefficients Ai and A2 appearing in ( |4.24p by 
considering homogeneous configurations as we have done here, because the corresponding 
terms involving derivatives of the velocity like a^ v and uj^ vanish, as the velocity remains 
constant. We can determine these coefficients by considering inhomogeneous configurations 
and adapting our method, provided we can argue that the metric should be regular in the 
ingoing Eddington-Finkelstein coordinates. 



5 Outlook for RHIC and ALICE 

We briefly mention here the relevance of this work for the great challenge of modelling the 
space-time evolution of the matter formed by ultra-relativistic heavy ion collisions at RHIC 
and ALICE. Experiments at RHIC suggest the validity of the following picture |38| : (i) a 
large fraction of the initial kinetic energy of the colliding ions is thermalized astonishingly fast 
(in time < 1 fm) forming a locally equilibrated hot and dense fireball parametrized by a profile 
of the hydrodynamic variables - namely the temperature, four- velocity and chemical potential 
fields, (ii) the strongly interacting fireball undergoes hydrodynamic expansion which can be 
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described by the hydrodynamic equations like the Navier-Stokes' equation, and (iii) the initial 
transverse hydrodynamic flow at the time of local thermalization in most cases vanishes. 
Most of the data at RHIC is in good agreement with this simplistic picture especially in the 
mid-rapidity region, i.e. for the most central collisions at the highest beam energy of y / Sjv r /v = 
200 GeV. However, despite the success in explaining the transverse momentum spectra of 
hadrons, the elliptic flow coefficient, etc., this does not reproduce pion interferometric data 
like HBT radii leading to the well-known RHIC HBT puzzle. 

It is necessary to have a better phenomenological model for the space-time evolution 
of the fireball to explain the data completely, and also to reduce theoretical uncertainties. 
The best theoretical tool at hand for studying evolution of strongly coupled matter of gauge 
theories in real time is the AdS/CFT correspondence. It is well-known that the AdS/CFT 
correspondence gives r]/s = l/4ir ||, while the current analysis of experimental data suggests 
1 < 4:7r\r)/s\ < 2.5 for temperatures probed at RHIC PU| . 

Here we will propose that the AdS / CFT correspondence can be used to develop a com- 
plete phenomenology for the evolution of the strongly coupled matter, describing both the 
late stages of local thermalization and the subsequent hydrodynamic expansion in an unified 
framework. These phenomenological equations involve a closed set of equations for evo- 
lution of the energy-momentum tensor and the baryon number charge current alone. For 
simplicity, if we assume that the baryon number chemical potential is zero throughout the 
evolution, then the full set of phenomenological equations are QL21| ), ( |4.22| ), ( |4.24| ). It should 



be straightforward to include the dynamics of the charge current but we will not attempt 
this here. 

The advantage of our proposal is that there is a very natural way to connect the expansion 
of the fireball with any model which describes the early stages of the collision process, as for 
instance the parton cascade model [EDI. All that is needed is to match the evolution of the 



energy-momentum tensor and conserved charge currents before and after the matter enters 
in the strongly coupled phase of evolution. Importantly, the matching with the initial regime 
does not require the energy-momentum tensor to be hydrodynamic. 

To state more concretely, the space-time evolution of the matter till chemical and kinetic 
freezeout with subsequent hadronization can be divided into two phases : 

• The entry into strongly coupled phase : This phase describes the very initial stage 
of particle production at high energies and can probably be described in the perturba- 



tive framework as in the parton cascade model |l0j . The energy-momentum tensor can 
be computed by such models and the local temperature fields can be obtained from the 
Landau-Lifshitz decomposition of the energy-momentum tensor. As mentioned before, 
any arbitrary energy-momentum tensor can be written in the Landau-Lifshitz form 
and thus we can always obtain the local effective temperature. However, this requires 
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using an appropriate equation of state for the initial phase. The effective temperature 
will automatically give us when the matter locally enters the strongly coupled phase. 
A more detailed analysis of this initial phase is beyond the scope of this paper. 

The fireball formation and its expansion : Once the matter enters into the 
strongly coupled phase, we can use the phenomenological equations from AdS/CFT to 
describe both the local thermalization and late hydrodynamic expansion in an unified 
framework. The initial conditions can be obtained by matching the energy-momentum 
tensor with the model describing the initial stage Q As the matter is expanding and 
boost invariance is a good approximation for central collisions, we may use the boost 
invariant versions of eqs. ( |4.21fl , ( f4.22p , ( |4.24| ) . 



Homogeneous relaxation in the boost-invariant regime which estimates the time for 
local thermalization or fireball formation has been studied in the linearized approxima- 
tion in [4~I|. Furthermore, as also shown in |5T], the homogeneous quasinormal modes 



in the static case can be mapped readily to boost invariant geometries. The boost- 
invariant version of our equations, as discussed here, must reproduce these linearized 
fluctuations of the boost-invariant geometry which are homogeneous in transverse co- 
ordinates. The lowest static homogeneous quasinormal mode gives an estimate of the 
time for local equilibration, i.e. fireball formation. In the boost invariant case, the 
time-dependence in the linearized plane wave approximation is not exponential, but 
given by a Bessel function as in y/rJ ± 3 (^uri^nTf)^^ Y where u ~ 2.74667(7rT/) and 
Tf is the fireball temperature which is about 175 MeV at RHIC. At late time, this 
reduces to a proper-time damping of the form : exp(^ — | • 2.74667 ■ (nTfr)% ■ ^ . 

In future work, we would like to study the boost-invariant version of the non-linear 
homogeneous equations given by ( |4.24| ) and see numerically how this behavior both 
qualitatively and quantitatively changes due to non-linearities. Further, it should be 
also possible to include inhomogeneities in transverse directions and also study the 
transition of the fireball to the hydrodynamic regime using the full set of equations 
- ( }4.21| ), ( }4.22| ), ( f4.24j ). In practice, one needs to use these equations in a coordinate 
system (as in 0) better adapted for the late equilibrium state of the fireball which 
is an ideal fluid undergoing boost-invariant expansion. This coordinate system com- 
prises of the proper time coordinate r of the late time expansion, the coordinate y 
parameterizing rapidity, and the two transverse coordinates x l and x 2 . 



12 This matching is however complicated as the equations of evolution involve time derivatives of arbitrary 
orders. 
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6 Discussion 



From the point of view of gauge/gravity duality, our work leads a few open questions. To 
conclude, we will discuss two of the most immediate and important of them here. 

• Is there an entropy current for homogeneous relaxation? In the case of pure 
hydrodynamic behavior, it is known that the horizon in the bulk metric also gets 
deformed mildly, and the deformation can be calculated order by order in the derivative 
expansion. This has been shown in [fl2| , where it has been also demonstrated that the 
pull-back of the Hodge-dual of the area form on the horizon world-volume to the 
boundary gives the construction of a family of entropy currents whose divergences are 
positive-definite, order by order in the derivative expansion. 

This is an interesting finding, because to our knowledge, this gives the first construction 
of an exact entropy current in the hydrodynamic regime, beyond any approximation 
like the limit of validity of the Boltzmann equation. An interesting question is whether 
gauge/gravity duality implies the existence of such entropy currents at large N, even 
beyond the hydrodynamic regime. 

In our case of spatially homogeneous relaxation, even with the explicit metrics obtained 
here, this is a question which is not easy to answer. This is because the location of the 
horizon r(t) as a function of time will fluctuate only perturbatively in the amplitude 
expansion, but its time derivatives will not be under control. Since we have summed 
the time derivatives only recursively here, this does not reveal the global nature of the 
horizon even perturbatively in the amplitude expansion. We would have to tackle this 
question possibly numerically later. 

Alternatively, we can also analyze the general phenomeno logical eqs. ( }4.21|) , ( |4.22| ) 



and ( |4.24|) and find conditions under which these imply existence of entropy current (s). 
Then, one can verify if these conditions are satisfied by the phenomenological coeffi- 
cients obtained by gauge/gravity duality. 

What are the field-theoretic definitions of the general phenomenological 
coefficients? Recently, it has been possible to obtain a systematic procedure for 
constructing Kubo-like formulae for all hydrodynamic transport coefficients (including 
those which appear only non-linearly) in terms of the low frequency and large wave- 
length expansion of n-point correlation functions of the energy- momentum tensor |[43| . 
It has so far not been possible to define such general formulae for non-hydrodynamic 
parameters beyond any approximation scheme. 

Our phenomenological coefficients (including those involved in homogeneous relax- 
ation), can be in principle obtained by constructing conservative solutions of non- 
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equilibrium field theory. However, we have been able to construct these solutions in 
fl5|1 in the perturbative regime in non-Abelian gauge theories only. It would be in- 



teresting to see if we can also find a non-perturbative method of constructing these 
conservative solutions and relate the phenomenological coefficients to n-point correla- 
tion functions of the energy-momentum tensor. 

If the above is possible, it would also provide a non-trivial consistency check of gauge/gravity 
duality. The n-point correlation functions of the energy-momentum tensor can be inde- 
pendently calculated by gauge/gravity duality, so we can obtain the phenomenological 
coefficients through them and see if they match with the values required by the regu- 
larity condition on the future horizon. 

• Is the ingoing Eddington-Finkelstein coordinate system the right choice for 
all configurations which equilibrate? It would be certainly of interest to see if 
the bulk metric is manifestly regular at the future horizon in the ingoing Eddington- 
Finkelstein coordinate system, for general configurations which equilibrate at late time. 
If this is so, we can readily extend our proof to show that the general regularity 
condition for the future horizon in the bulk, is indeed given by ( 4.21 ), ([4.22 ) and 
( f4.24| ), for the right values of the phenomenological coefficients. 
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A The Feffer man- Graham vector and scalar equations 
at the second order in the amplitude expansion 

The vector eq. ( j2.6| ) at the second order in the amplitude expansion reduces to 
tufv^ + Z-^ f^^u,^) = uV lr (p)fW(p,co 1 )f^(p,u-uj 1 ) 

+V 2r (p)(uj + wi)/«0o, u - ^i)^/ (1) (P ; "0 
+V 2r (p){2uj - Wl )/«(p, u; x )|VW - wi),(A.l) 
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where 

12 P 3 T , , \ 4p 3 6 4 (46 4 + 3p 4 ) , « 6 4 , A , 

Vl{p) = W^7~Y Vlrip) = W^7)WT7T 2 > VM = W^7)- { ] 

A close inspection shows that the right hand side of ( |A.1|) is symmetric under the exchange 
of u\ with oj — L)\ as it should be. 

The scalar eq. ( |2.7| ) at the second order in the amplitude expansion reduces to 

d d 2 \ 

S u (p) + S Mq^ + Ssl ^Q^ J /a (P> w > w i) 

(d d 2 \ 

= S lr (p)f W (p, ^i)/ (1) (P^ - wi) + S 2r (p)^ 1 \p,co 1 )^ 1 \p,co - Ul ) 
+SM (p, ui) ( P ,cj-u 1 ) + /« (p, w - w,) A/W (p, Wi) 
+^ r (p) (/ (1 Hp,^i)^/ (1) (p^-u; 1 ) + /( 1 )(p,a;-u; 1 )^/ (1) (P^i) ), (A.3) 



where, 



_ . \ 24p 2 (46 4 -p 4 ) 8p 2 (1926 12 + 48p 4 6 8 + 36p 8 6 4 + r 12 ) 

^lKPj = 7771 , ITT - > ^4«(P) " 



(4M + p 4 ) 3 ' " vry (46 4 -p 4 ) 4 (46 4 + p 

3(46 4 + 5p 4 ) . . (-166 8 + 48p 4 6 4 + 5p 8 ) 

MP) = -77U I „dV> » 6 5/(p) - 



4^1 



p(46 4 + p 4 ) 2 ' utvrv p(46 4 -p 4 ; 



3 



O f N 3 N (4& 4 + p 4 ) 

^(P) = -pi^, ^(P)= (4&4 _ p4)2 - (A.4) 
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